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PRECISION DETERMINATION OF NUCLEAR GYROMAGNETIC 
RATIOS! 


By E. W. GupTILL, W. J. ARCHIBALD, AND E. S. WARREN 


Abstract 


The ratio of the Larmour frequency of precession for hydrogen and fluorine, 
and for sodium and aluminum have been measured with an accuracy of one part 
in 105: 

Yi /’p = 1.062917 + 0.00001 
Yna/¥a1 = 1.015081 + 0.00001 


These ratios are corrected for the atomic diamagnetism giving for, y, the 
nuclear gyromagnetic ratios: 


Yq/Vp = 1.06245 + 0.00002 
Yya/Va) = 1.01494 + 0.00002 


A brief description is given of the relevant equipment used in obtaining these 
measurements. 


Several methods for measuring nuclear gyromagnetic ratios have been 
developed in the past few years (1, 7, 8). The present article describes a 
slightly different method which lends itself to an accurate determination of the 
ratios of nuclear moments. The technique involves the use of a single sample 
in which both nuclei are to be found. The measurement of a ratio has been 
reduced to the measurement of the Larmour frequencies of the two nuclei, 
both resonances being detected by the same electronic device, the change from 
one resonance to another being effected with no physical adjustments other 
than the resetting of an oscillator. 


Apparatus 


Fig. 1 contains a block diagram of the apparatus. Element A of the diagram 
is a combination oscillator and detector, a detailed circuit diagram of which 
is to be found in Fig. 2. The coils which in Bloch’s experiment are called the 
driver coil and pickup coil are in this apparatus respectively the plate coil 
and grid coil of the oscillator. The amplitude of oscillation of A is determined 
(among other things) by the mutual inductance of these two coils, which is 
altered by the nuclear precession when it occurs. 


The magnetic field is modulated about 0.0025% by the 30 cycle oscillator D. 
This slight modulation, as in Purcell’s method (2), is a small fraction of the 


1 Manuscript received January 31, 1950. 
Contribution from the Department of Physics, Dalhousie University, Halifax, N.S. 
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natural line width. When in the region of the resonance, there is a 30 cycle 
amplitude modulation of the oscillator A. The oscillations of A are detected, 
giving a 30 cycle signal which is amplified by B, which is in turn detected by 
the phase-sensitive detector (3) C. The 30 cycle oscillator D is coupled to C 
in order that only those signals having the same frequency as D may give rise 
to a unidirectional current through the meter E. 
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Fic. 1. Block diagram of apparatus. 


It is essential that the perturbing radio-frequency field at the nuclei be 
small if a net nuclear magnetic moment is to be preserved. The control of this 
r.f. field can be achieved in two ways: 

(a) by altering the mutual inductance between the grid and plate coils; and 
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(b) by varying the grid bias of the semi-remote cutoff tube 6SG7 


(Cf. Fig. 2). A point by point plot of the meter reading E against the mean 
magnetic field yields either the slope of the dispersion or the absorption curve 
or a superposition of both, depending upon the phase change suffered by the 
signal induced in the grid coil. In the apparatus as herein described the phase 
was adjusted by JL, in Fig. 2 to give the slope of the dispersion curve. The 
shape of the line may also be altered (but to a lesser extent) by the resistor in 
the cathode circuit of the 6SC7. , 


Magnet 


The magnetic field was provided by a one-ton magnet, the iron for which 
was presented to the University by the Steel Company of Canada. The ratio 
of pole face diameter (8 in.) to gap width (1 in.) was very close to that of the 
McGill cyclotron magnet, and the shimming was copied directly from it and 
machined into removable pole faces. It should be noted that the field require- 
ments for our magnet are, strictly, not the same as for a cyclotron, the former 
requiring a field which is as homogeneous as is possible. However, the shimming 
adopted was satisfactory for fields below 5000 gauss. Fig. 4 shows a plot of the 
percentage change in the field versus radial distance from the center. The 
data for this curve was taken at 3800 gauss. In the present experiment only 
extremely small drifts in the magnetic field could be tolerated. It was, con- 
sequently, necessary to control the current electronically with the circuit 
shown in Fig. 3. This particular control is intended for a constant voltage 
source such as storage cells, and with it the drift in magnetic field could be 
kept as low as 2 or 3 parts in 10° over the period of an hour. The resistance 
of the magnet wire would change about 1% in a corresponding period. 


Experimental Results 
Hydrogen and Fluorine 
A sample containing about 0.5 cc. of hydrofluoric acid was used for the 

hydrogen and fluorine measurements. The variation in magnetic field over 
the volume of the sample was about 5 parts in 10°. It was assumed that this 
was comparable to the natural line width in this sample, since the actual line 
width was reduced very little by a further reduction of the volume. The signals 
from both the hydrogen and fluorine were large enough to be viewed directly 
on the screen of the oscilloscope when a 60 cycle modulation was applied to 
the magnetic field. The procedure was to adjust the frequency of oscillator A 
until the signals from the hydrogen nuclei were in the center of the oscilloscope 
trace. The frequency of the oscillator was then measured. A similar measure- 
ment was made for the fluorine nucleus, the only alteration in the equipment 
being the condensers controlling the frequency. Twenty-one measurements of 
vy/vp Were obtained—12 at a field of 3500 gauss and 9 at 1000 gauss. 

3500 gauss: 1.062916 + 0.00001 

1000 gauss: 1.062919 + 0.00001 

Average: 1.062917 + 0.00001. 
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The second group of measurements was taken after altering as many variables 
as was possible, including the oscillator A and all its components. This was 
done in an effort to assess the magnitude of systematic errors. The estimated 
error (0.00001) is 1/5 of the line width and 4 times larger than that obtained by 


; 2d? ; : oi 
evaluating Ne where JN is the number of ratios and d the deviation from 


the mean. Seventy-five per cent of the ratios fell within the limits of the given 


error. 
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Fic. 4. Graph showing the variation in mignetic field near the center of the magnet. 


Sodium and Aluminum 


One-half a cubic centimeter of sodium aluminate was used for the measure- 
ments on Na and Al. As before, the line width was comparable to the field 
inhomogeneity; however, the signal-to-noise ratio was smaller, and it was 
necessary to use the 30 cycle modulation and the lock-in amplifier. Fig. 5 shows 
an Esterline—Angus plot of these lines as the magnetic field was slowly varied 
through resonance. A typical ratio was measured by first setting the oscillator 
frequency equal to the Larmour frequency of the magnetic moment. The 
criterion for the frequency adjustment was the maximum deflection of the 
Esterline-Angus meter. The frequency was then measured and similar data 
taken for the second magnetic moment. Thirty values of this ratio were 
obtained at a field of 3400 gauss: 


Yna/VAI = 1.015081 = 0.00001 


In this case the distribution was not gaussian but rather with one maximum 
each side of a central maximum in which 40% of the values occurred. The 
reason for this is not clearly understood, but there are good indications that 
it is a function of the apparatus. The mean of the 12 values occurring in the 
central region is 2 parts in 10° lower than the ratio given which is the mean of 
all values. The spread of the two subsiduary maxima was 3 parts in 105. 
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The frequency ratios above may be corrected to yield Y¥y/Y¥r and Yy4/Yai, 
where the y’s are the nuclear gyromagnetic ratios. This diamagnetic correction 
has been calculated by Lamb (6): 


Hy = Hy) —— | in) g 
3mc Jr r 
where H, is the diamagnetic field at the nucleus, Hp the magnetic field in which 


the charge density f(r) is placed, e and m the charge and mass of the electron. 
The integral can be evaluated for any particular atom by the Hartree—Foch 
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Fic. 5. Esterline-Angus (recording milliammeter) plot of signals from aluminum and sodium 
in sodium aluminate. There is a discontinuity (1.5%) in the magnetic field at the center of the 


trace. 


method of the self-consistent field. In the case of the hydrogen atom H;/H)= 
1/3a?= 0.177 X 107‘, where a is the fine structure constant. For the hydrogen 
in H,O the factor would be almost twice this, or about three parts in 10°. 


The following values of H;/H») were computed from the wave functions 
given by Kennard (Na) (5) and Hartree (Al and F) (4); Na 0.63 X 107%; Al 
0.77 < 107%; and F 0.47 X 107%. These corrections, give for the nuclear 
gyromagnetic ratios, 

Yu/VYr = 1.06245 + 0.00002 
Ynwe/Vai = 1.01494 + 0.00002. 
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SPECIFIC HEATS OF CERTAIN SALTS OF IRON GROUP 
ELEMENTS FROM 65° TO 300°K.! 


By H. D. VAsILerF? AND H. GRAYSON-SMITH?® 


Abstract 


Using a new low temperature calorimeter, which is briefly described in the 
paper, the specific heats have been measured from.65° to 300°K. for the following 
salts: chromium sulphate (hydrated and anhydrous), chromium nitrate, cobalt 
nitrate, and nickel nitrate (hydrated). Hydrated chromium sulphate was found 
to have a transition of the second order at 195°K., while the specific heat of the 
anhydrous salt was quite regular. The hydrated nitrates all showed second 
order transitions in the neighborhood of 150°K. The entropy changes associated 
with these transitions have been estimated approximately, and vary from about 
0.4R for cobalt nitrate to 1.65R for chromium nitrate, where R is the gas con- 
stant. Pending further evidence, it is tentatively suggested that the transitions 
are due to the onset of partial rotation of the H:O groups in the crystals. 


Introduction 


As part of a study of simple compounds of elements of the iron group, 
chromium to nickel, the specific heats have been measured from 65° to 300°K. 
for the salts listed in Table I. The same salts have been investigated magneti- 
cally by Johnson (8), using materials from the same lots. The two investigations 
were carried out at the same time because knowledge of either property may 
help in the interpretation of the other. However, the magnetic properties 
proved to be quite regular over the temperature range concerned, while the 
calorimetric study has revealed that all the hydrated salts pass through tran- 
sitions of the second order. These are discussed in the present paper. They 
must be due to some cause other than magnetic, but cannot be definitely 
identified on the calorimetric evidence alone. 


Materials Used 


The salts were obtained from British Drug Houses (Canada) Ltd., and were 
specified as having impurity less than the amounts listed in Table I. The 
anhydrous chromium sulphate was prepared from the hydrate by heating at 
400°C. until the weight was constant. The observed loss of weight was 1 to 2% 
greater than the amount accounted for by the removal of 18 H,O. In a salt 
such as Cre(SO,4)3.18 H,O there is always some uncertainty concerning the 
amount of water of crystallization actually present. From the loss of weight 
in heating it was concluded that the amount of water in the sample used was 
close to the theoretical value. 


1 Manuscript received February 7, 1950. 

Contribution from the Department of Physics, University of Toronto, Toronto, Ont. 
2 Present address: Defense Research Board, Oltawa, Ont. 
3 Present address: Department of Physics, University of Alberta, Edmonton, Alberta. 
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TABLE I 
SPECIFICATIONS OF SALTS INVESTIGATED 








Chromium sulphate, Cr2(SO4)3.18H20 Cobalt nitrate, Co(NOs3)2.6H20 
hloride (Cl)...... 0.003% Chloride (Cl)...... 0.001% 
PO APO) 56s so Kce oe 0.01% Sulphate (SOx,)..... 0.01% 
ike tree 0.5% Po eee 0.003% 
Nicke! (Ni). ......5..- 0.002% 
Chromium sulphate anhydrous, Cr2(SOu)s Alkalis and alkaline 
Same as hydrated salt earths (Na)......0.03% 
Ammonia..... No reaction. 
Chromium nitrate, Cr(NO3)3.9H2O 
Chloride (Cl)...... 0.003% Nickel nitrate, Ni(NO)3)2.6H20 
Sulphate (SO4)..... 0.01% Chloride (Cl)...... 0.002% 
PEOHRO).3 5. 5 eae 0.01% Sulphate (SO,)..... 0.01% 
PRIIB S56) Ssn coe 0.5% Cobalt (CO)... 25. 0.0005% 
BOOETEO No.0 coe con 0.001% 
Alkalis and alkaline 
earths (Na)......0.038% 





The Calorimeter 


The low temperature calorimeter is shown in Fig. 1. It was of similar design 
to that used by Elson, Grayson-Smith, and Wilhelm (2), but was completely 
rebuilt. The powdered samples were compressed into the calorimeter chamber 
E by means of the plate G and nut H. The cover plate J was then soldered 
into place with Wood’s metal, the calorimeter was exhausted, and filled with 
helium gas at a few centimeters pressure in order to improve the thermal 





j 
| 
| 
| 





Fic. 1. Low temperature calorimeter. 
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contact at low temperatures. The calorimeter was heated from the outside 
surface by means of a constantan coil K, protected by a radiation shield D. 


The temperatures were measured by means of resistance thermometers A, B, 
wound on a mica frame, and mounted in the central tube F, which was also 
filled with helium. The principal thermometer was of platinum, with a resis- 
tance at 0°C. of 8.4782 ohms. For temperatures down to 91°K. it was cali- 
brated according to the specifications for the International temperature scale. 
At lower temperatures Henning’s deviation method (6, 7) was used. Using a 
Mueller bridge, and a galvanometer with a sensitivity of 0.2 wv. per mm., 
resistances could be measured to 10~‘ ohm, and temperature changes to 0.003°. 
The increases of temperature in the specific heat determinations were 0.5 to 1°, 
and could therefore be measured with a precision of 0.6 to 0.3%. 


Auxiliary manganin thermometers were wound on the surface of the calori- 
meter at C, along with the heating coil, and on the inner mica frame, along 
with the platinum thermometer. The inner manganin thermometer was intended 
for extremely low temperatures, where the platinum thermometer becomes in- 
sensitive. It was not used in the measurements to be reported in this paper. 
That on the outer surface was used in tests of the rate at which temperature 
equilibrium was attained in the calorimeter. 


The rate of heating was determined by measuring the potential difference 
across the constantan heater with a Leedsand Northrup type K2 potentiometer. 
The resistance of the heater was determined by comparing this potential drop 
with that across a standard resistance connected in series in the same circuit. 
This needed only to be checked occasionally, since the resistance of the con- 
stantan was nearly constant over a 1° temperature interval, and could be taken 
equal to the value at the middle of the interval. The time of heating, and the 
times at which temperature measurements were made, were recorded on a 
tape chronograph, and compared with time marks made by a 25 cycle electric 


clock motor. 


The whole calorimeter assembly shown in Fig. 1 was suspended by threads 
inside a protecting jacket, which could be evacuated for insulation, or filled 
with a gas for thermal contact. Whenever possible, the jacket was kept at a 
constant temperature by immersing it in a cryostat of liquid nitrogen, liquid 
oxygen or liquid ethylene. For temperatures between 91° and 135°K., between 
the oxygen and ethylene ranges, the mean temperature of the jacket could 
be held sufficiently steady to enable measurements to be made by immersing 
it to different depths in liquid oxygen. For temperatures above the boiling 
point of ethylene (170°K.) the calorimeter jacket was suspended in the stream 
of cold air rising from liquid air boiling at a steady rate. 


In low temperature calorimetry the influx of heat from parts of the apparatus 
at room temperature, by conduction along the electrical leads, and by radiation 
through the pumping tube, is usually greater than the loss of heat by radiation 
to the jacket, which is at the cryostat temperature. This source of error was 
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reduced by wrapping all the electrical leads around copper posts attached to 
the calorimeter jacket, and by placing a polished copper baffle over the outlet 
to the pumping tube. 


Correction for Heat Exchange 


When low temperature specific heat measurements have not been made in 
an adiabatic calorimeter it has been usual to correct for heat exchange with the 
surroundings by the graphical method developed by Keesom and Kok (9). 
This method involves three approximations, which may or may not be satis- 
factory for a particular calorimeter. It is assumed that the loss of heat to the 
immediate surroundings is proportional to the temperature difference, that 
the rate of gain from the more distant surroundings at room temperature is 
constant, and that the temperature of the calorimeter and sample is uniform. 
In preliminary tests the form of the temperature-time curve before and after 
a heating interval showed that the first two approximations were sufficiently 
accurate to be used in correction for heat exchange. However, during and 
immediately after a heating interval there was an appreciable difference of 
temperature between the center and the surface, such that, with heating 
intervals of two to four minutes, the maximum temperature at the center was 
not attained until 25 to 50 sec. after the surface heater was turned off. 


In order to determine whether this time lag would cause an appreciable 
error in the temperature correction, a theoretical calculation was made, based 
on the well known equations for a cylinder with a radial temperature gradient. 
The conclusion was that the graphical method of Keesom and Kok would give 
too high values of the heat capacity when the time lag was appreciable. Cor- 
rection for the time lag was not feasible on account of the uncertainty in the 
mean thermal conductivity within the calorimeter. However, any specific heat 
measurement involves the difference between the heat capacities with and 
without the sample. If these are both measured experimentally by the same 
procedure, the errors due to time lag should approximately cancel, and the 
specific heat values should not be seriously affected. 


Heat Capacity of the Empty Calorimeter 


The heat capacity of the calorimeter alone was measured directly, and was 
also deduced from a series of measurements with a copper sample, using the 
values given by Dockerty (1) for the specific heat of copper. The results are 
shown in Fig. 2, where the circles represent the direct values, and the triangles 
indirect values. The agreement between the two sets of values shows that the 
time lag has not seriously affected the differences, in spite of the large difference 
in thermal conductivity between copper and helium gas. 


From the scattering of the points in Fig. 2 it was estimated that the random 
error in subsequent measurements of the specific heats of the salts studied 
should not exceed 0.001 cal.-gm.~!-deg.—! 
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Fic. 2. Heat capacity of calorimeter. ©, direct measurements; A, deduced from measure- 
ments with copper sample. 


Specific Heats of Iron Group Salts 


Measurements of the specific heats of the salts listed in Table I are shown in 
Figs. 3 to 6. In each case the plotted points represent at least two distinct 
series of measurements. Values obtained from the smoothed graphs are col- 
lected in Table II. 


The most interesting feature of these results is that all the hydrated salts 
show specific heat anomalies, while the specific heat of the anhydrous chromium 
sulphate is quite regular over the range of temperature investigated. From 
the shape of the specific heat curves it is evident that the anomalies are due 


to transitions of the second order. 


’ The increase of entropy from 65° to 273°K. was determined by drawing 
smooth curves of the ratio (heat capacity per mole)/(absolute temperature) 
and integrating this function numerically, with the results given in Table III. 
It is not possible to extrapolate the entropies to 0°K. because these salts are all 
paramagnetic, and may be expected to have magnetic transitions of some kind 
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TABLE II 
SPECIFIC HEATS DERIVED FROM SMOOTHED GRAPHS 


Specific heat, cal.-gm.—!-deg.— 
Temperature, 
°K. 

















Chromium Chromium . : 
Chromium Cobalt Nickel 
ioe leans nitrate nitrate nitrate 
300 0.312 0.172 0.273 0.371 0.382 
280 0.304 0.163 0.255 0.363 0.365 
260 0.293 0.154 0.237 0.350 0.345 
240 0.280 0.146 0.216 0.333 0.323 
220 0.265 0.137 0.193 0.311 0.297 
200 0.252 0.128 0.169 0.283 0.269 
195 0.312 - - - - 
190 0.273 - - - - 
185 0.251 - - = -~ 
180 0.239 0.120 0.151 0.250 0.235 
170 0.224 - 0.149 0.234 0.217 
160 0.214 0.111 0.168 0.240 0.212 
155 - - 0.182 0.270 0.222 
150 - = 0.167 0.275 0.244 
145 - ~ 0.156 0.256 0.237 
140 0.201 0.102 0.150 0.248 0.220 
130 - - 0.143 0.236 0.201 
120 0.181 0.092 0.139 0.224 0.188 
100 0.151 0.080 0.133 0.192 0.160 
90 0.136 0.072 0.129 0.172 0.144 
80 0.120 0.064 0.122 0.150 0.128 
70 0.104 0.053 0.110 0.125 0.110 
65 0.095 0.046 0.100 0.112 0.102 | 
| 
TABLE III 
ENTROPY CHANGES AND TRANSITION TEMPERATURES 
Total entropy change Transition Estimated entropy 
Salt Sor3— Ses temperature, due to the transition 
cal.-mole.—! -deg.—! oe (as multiples of R) 
Cr2(SO4)3.18H,0 202.5 195 1.41R 
Cr2(SO,4)s 56.2 — a 
Cr(NO3)3.9H:0 88.2 158 1.65R 
Co(NO3)2.6H:O 93.5 153 0.4R 


Ni(NO3)2.6H2O0 87.1 149 0.6R 


at lower temperatures. Table III gives also the transition temperatures, and 
the estimated increases of entropy due to the transition (as multiples of the 
gas constant R). For the cobalt and nickel nitrates the entropy change due 
to the transition can only be roughly estimated because the specific heat of 
the crystal lattice is also changing rapidly over the temperature range con- 
cerned, and seems to have a discontinuity at the transition temperature. 


Evidence was obtained that supercooling effects might accompany any of 
these transitions. In most of the experiments it took several hours to cool the 
calorimeter and sample to the required initial temperature for a series of 
readings, using gas conduction only. Under these conditions the specific heat 
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curves in the neighborhood of the transitions were reproducible. However, if 
the salts were quenched rapidly to 90°K. before a series of measurements, 
abnormally high values of specific heat were obtained between 90° and the 
transition temperature. This indicates metastable forms of higher entropy, 
and so more or less ‘‘frozen-in disorder.’”’ After chromium sulphate hydrate had 
been quenched, it had to be annealed at room temperature for several days 
before the original curve could be reproduced. 


Discussion 


Second order transitions of the type described are well known, and may be 
due to any one of a number of causes. In pure substances they have been 
ascribed in different cases to ferromagnetism, to ferroelectric effects, or to the 
onset of partial rotation of one or more molecular groups within the crystal. 
Additional evidence must therefore be obtained in order to identify the cause 
in the salts investigated. It is known that they are not magnetic, since Johnson 
(8) has measured the paramagnetic susceptibilities of the same materials, and 
has found them to be quite regular down to 77°K. The possibilities remain of 
either dielectric transitions similar to that observed by Roberts (12) and by 
Harwood et al. (5) in barium titanate, or of rotational transitions similar to 
those observed in a number of ammonium salts (4). 


Of these alternatives the onset of rotation seems the most likely. Since the 
anomaly occurs with hydrated chromium sulphate and not with the anhydrous 
salt, it must, in that case at least, be due to the water of crystallization. So 
far as the authors are aware, rotation within a crystal of H,O groups has not 
hitherto been suggested as a cause of second order transitions. The nearest 
parallel case in the literature is a transition in Ni(NO3)2.6NHs, which has been 
ascribed by Long and Toettcher (11) to the onset of rotation of both the NO; 
and NH; groups. Unexplained transitions have been observed in lanthanum 
magnesium nitrate (3) and hydrated magnesium chloride (10), and might be 
due to the same cause as those reported in the present paper. 


Whether or not this suggestion is correct cannot be established on the 
calorimetric evidence alone. The obvious next step is a study of the dielectric 
properties. Apparatus for this study is now being constructed. 
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MATHEMATICAL ANALYSIS OF AN ACOUSTIC FILTER! 


By N. OLson 


Abstract 


Equations have been derived for the frequency limits of the attenuation bands, 
the degree of attenuation, and the phase change of a specific type of filter which 
cannot be readily: analyzed by the side branch method. The characteristics 
have been verified by measurements on experimental filters. 


Introduction 


A theory of acoustic filters has been developed by Stewart (3, 4) analogous 
to the theory of electrical filters, in which the series and shunt impedances 
of electrical filters are replaced by appropriate.lumped acoustic impedances. 
In particular, a conduit with regularly spaced side branches may be repre- 
sented thus, and the formulas for electrical filters can be used. In this theory, 
the assumption is made that the wave length of the sound is much greater 
than the dimensions of the filter section, so that the phase change in each 
section may be neglected. 


A later theory developed by Mason (1, 2) considers the effect of the side 
branches on the transmission of a plane wave through a conduit. For the 
special case in which the wave length of the sound is much greater than the 
dimensions of the filter, this theory reduces to that of Stewart. Mason also 
extends his theory to include the case of a cylindrical side branch, coaxial with 
the main conduit and coupled to it via a gap in the conduit. The side branch 
impedance is determined by assuming that a radial piston acts in the opening. 


This paper deals with a specific type of filter consisting of a conduit with a 
series of expansions at regular intervals. By considering the transmission of a 
plane wave through the conduit, equations will be derived for the limits of the 
attenuation bands, phase change, and attenuation. The assumption is here 
made that the transmitted wave is plane at all points along the filter, an 
assumption which is valid if the change in cross section is not excessively large. 
This type of filter can be considered as a special case of the one discussed by 
Mason, with the side branch opening extended until it includes the entire 
length of the side branch. However, Mason’s approximations are not appli- 
cable to this case. 


Filter Theory 


The filter is taken to consist of an infinite series of sections as shown in 
Fig. 1, (a), or 1, (0), and in which the impedance is the same at all junctions 
1 Manuscript received February 8, 1950. 
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Fic. 1. Filter sections. (a) m = Si/S2< 1. (6) m = Si/S2> 1. 


between sections. 5S; and S» are cross-sectional areas. The velocity potential 
in a plane wave is given by 


¢ as Ae’ t— + Bei + kx) (1) 
from which the excess pressure and volume current are 

p = — ipokc(Ae~*** + Be*k*)e' 

X = —ikS(Ae~**? — Be*#*)ei, 


w = 2m X frequency, 

k w/c, 

c¢ = velocity of sound, 

po = equilibrium density of medium, 

S = cross sectional area, 

A and B are arbitrary constants, 

X= Sé, where — = velocity of particles in medium. 


Omitting the exponential in ¢, and starting at position 1 in the filter section, 
Fig. 1: 


(4) 


pr = seeds HO 
Xi = —tkS;(Ai — Bi)» 


At 2, 
p2 = —ipokc(Aye7 "2 + Bye**'/2) \ (5) 
X. = —ikSi(Are~#/2 — Byettl/2y J’ 
pb, = —ipokc(Ase "2 + eee? (6) 
X; = —ikS2(Ace~*’2 — Byet#l/2) 


where unprimed quantities represent conditions to the left, and primed quan- 
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tities represent conditions to the right of the junction between the two cross 
sections. 


At 3, 
Ps = —ipokc(Age*!/2 4. Ee (7) 
¥. 2 —ikS2(Agei*¥/2 Zee Bye**l/2) 
b; = —ipokc(Ase~**/24 Byeittl/ ot. (8) 
= 2 —ikS;(Ase#*Y2 — Bye**"/2) 

At 4, 
Pa = —ipokc(Ase~*! + Bye?) \. (9) 
Xs = —ikSy(Ase~**! — Bye**!) 


Since there is continuity of pressure and volume current at the junctions, 
Equations (4) to (9) may be reduced to 
fs = bi/4m|{ (m + 1)? cos 2k] — (m — 1)*} 
—iZ1/Zo{(m + 1)? sin 2k/ + 2(m? — 1) sin ki}] 
X,= X1/4m{{ (m + 1)* cos 2kl —(m — 1)?} . (10) 
—iZo/Z1{ (m + 1) sin 2kl — 2(m?— 1) sin k/}] | 


where p4/X4 = p1/X1 = Zo . Zo is “characteristic impedance” of the filter. 





Zi = pot/Si 
m= Si/ Se. 
From Equation (10), 
Zo~ fig/ ee Dee Ge ee 
(m + 1) cos kl — (m — 1) 


Eliminating Zo/Z, from Equation (10), 
1 2 al aa 1 2 
on. p,| + 1)? cos 2kl —(m ) 


4m 





~ / (m + 1)* sin?2kl — 4(m* — 1)? sin?kl 
16m (12) 
4m 
16m? 


X. = x, 


The two terms in each of the above expressions bear a simple relation to each 
other. If the first term is denoted by cos Y, then the second term will be sin Y. 


Thus, 
sin Y = sf Se sin2kl cos kl — (==) ’ (13) 
4m? m+ 1 


A ye SR 


| 
| 
' 
| 
| 
| 
| 
i 
| 
| 
| 
| 
| 
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ps = Pi(cos Y —7sin r) = ae (14) 
Xs = Xi(cos Y —isin Y)= Xie~*¥ 5 
An examination of Equation (13) shows that sin Y, hence Y, may be real or 


imaginary. Let Y = a + ib, where a and b are real, and represent phase angle 
and attenuation factor respectively. Then 


sin Y = sin a cosh b +i cosa sinh b. (15) 


Equations (13) to (15) determine the characteristics of the filter. Consider 
briefly the following three cases: 


(1) Cos? kl > (= ea iy. 
m+ 1 


Sin Y is real, and Y is real. Hence 6 = 0, and sin Y = sina. 


(2) Cos? kl< (2—)" 
m+ 1 


Sin Y is imaginary. Hence sin a = 0, and sin Y = 7 sinh b, from which Y =1b. 
Then from Equation (13) 


sinh b = fee sin? al (2) — cos? | ; (16) 
4m? m+1 





For attenuation, the negative value of 6 must be chosen. 


af 2 
(3) Cos? kl = (==) ,or sin Bl =.0. 
+1 


m 
Sin Y = 0, and sin a = 0, b = 0. 


Fig. 2 shows the manner in which the terms of the expression for sin Y vary 
with k/, hence frequency. The characteristics of the filter may be studied with 
reference to this, and Equation (14). As the frequency increases from zero to 


, cS m—1\? ‘ 
the first frequency for which cos? k/ = (2— , the attenuation factor re- 


m 
mains zero, and the phase angle a increases from 0 to z. This is the first trans- 
mission band. At frequencies for which + (==) > cos kl > — (=—1) ‘ 
m-+ 1 m+ 1 
there is attenuation, and a remains constant at 7. The attenuation factor 0 is 
determined from Equation (16). In the second transmission band, attenuation 
is again zero, and the phase angle a increases from 7 at the lower frequency 
limit, to 27 at the center, for which sin k/ = 0, to 3z at the high frequency limit 
of the band. 


It is seen that there will be an infinite series of successive transmission and 
attenuation bands. The band widths are determined by the value of m, and 














OLSON: MATHEMATICAL ANALYSIS OF AN ACOUSTIC FILTER 381 






a 2 
the cutoff frequencies by the expression cos’ kl = (==) . Maximum at- 
m 
| tenuation occurs at cos? k/ = 0. It will be noted that 
o=3 = (&=3) "2 where g = 1/m. 
m+1 g+1 






Hence the filters of Fig. 1, (a), and 1, (0), will have the same characteristics, 
except for the characteristic impedance Zo, Equation (11). 








Cos kl 





T, transmission bands; A, atten- 







Fic. 2. Graphical representation of filter characteristics. 
uation bands. 


Characteristic Impedance 






The value of Z) may be investigated in the same fashion as the transmission | 
characteristics. Rearranging terms in Equation (11), : 


Aa=21 / cos t+ (Fi) (17) 


V cos kl — ("= 1) | 






m+ 1 





m— 1 


m+1 
and therefore positive. Hence the positive root must be taken. In the atten- 
uation bands, Zo is a pure reactance, and may be positive or negative. At zero 


frequency, cos k/ = 1, and 
Zo= Vm Z,. (18) 





In the transmission bands, cos? k/> ( ) , Zo is real, a pure resistance, 






At cutoff frequencies, cos? kl = (2. Therefore, if m <1, Fig. 1, (a), 
m 
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—1\. : ; 
(= ) is negative, and Zp) = 0 at the first cutoff frequency, and increases to 
m 
infinity at the upper limit of the attenuation band. On the other hand, if 


m > 1, Fig. 1, (0), (* a ) is positive, and the zeros and infinities of Zo are 
m 


interchanged. Hence Zp can be traced through successive transmission and 
attenuation bands. The two cases are shown graphically in Fig. 3 for the actual 


ee 


ty 
SS) 
= 
c 
ql 
Q 
z= 
_ 
Y 
> 
° 
Y 
< 


a 
2 
a 
= 
° 
*) 
’ 
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po 


1000 00 5000 10000 
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3s 


Fic. 3. Impedance curves. T, transmission; A, attenuation. 
1. Characteristic impedance Zo, filter 6a, m = z 
. Characteristic impedance Zo, filter #6b, m > 
. Output coupling impedance Zao, filters #6a and ts, l= 3.38 cm., Z = poc/S=9.8. 
. Characteristic impedance Zo, filter #7a, m< 1 
. Characteristic impedance Zo, filter #7b, m> 1 
5. Output coupling impedance Zaza, filters #7a and #7, le= 6.7 cm., Z = poc/S=9.3. 


filters used in the measurements. The reactances are shown positive. In elec- 
trical filter theory, the section under examination may lie between successive 
“‘mid-series” or “‘mid-shunt”’ points in the infinite line on which the theory is 
based. The acoustic equivalents lie between successive ‘‘mid-contraction’”’ or 
‘“‘mid-expansion”’ points in the infinite acoustic line. In both systems, the fre- 
quencies of infinite impedance for one section are frequencies of zero imped- 
ance for the other, and vice versa. 


Impedance Matching 


The preceding theory applies to an infinite series of filter sections; i.e., each 
section is perfectly terminated. In the practical filter, with a finite number of 
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sections, the ends must be terminated in devices with impedances equal to Zo. 
As with electrical filters, exact impedance matching can be achieved only at 
certain frequencies. In electrical filters, the usual procedure is to terminate the 
filter in a pure resistance, equal to the characteristic impedance in the trans- 
mission band. In the present case, this could be achieved by terminating the 
filter in an infinite conduit, of cross section S, such that poc/S is equal to Zo 
in the transmission band. In practice, however, terminating elements of finite 
dimensions must be used. Since the application of the present development was 
simply to attenuate noise, a simple system was desirable. Hence, the effect 
of terminating the filter in short conduits was investigated. 


Za S 2, 


Fic. 4. Filter terminating conduit. 


Consider a short conduit of length /;, and cross section S, Fig. 4. The im- 
pedance Z, at the input end can be expressed in terms of Z, at the output end. 


Thus 


ee z| 2 iZ tan a ‘ (19) 


Z+1Z, tan kl; 


where Z = poc/S. Note that if Z,= Z, then Z,= Z also. The sections com- 
prising the filter are inserted between two conduits, Fig. 5, of lengths /; and /, 
and of cross section S, such that Z = Zo approximately, for the first trans- 
mission band. If the output conduit opens into free space, the impedance 
will be 


Rk , tpow 
Den ut See she ee 20 
b on + C, (20) 


where Cp, is the conductivity of the orifice. The first term represents radiation 
resistance, and will be small for wave lengths large compared to the dimensions 


OUT PU 


Zar S 


ds 


Fic. 5. Single filter section with terminating conduits. 
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of the filter sections. The second term is the inertance of the air in and around 
the opening, and at low frequencies this will also be small compared to Z. 
Therefore, as a first approximation, let Z,.= 0. Then from Equation (19), 


Z.2= iZ tan kl. (21) 


If now Za= Zo, then Lyi= 20= 1Z tan Rls. 
Substituting Z,; into Equation (19), 


Za= 1Z tan k(1i+ 2) . (22) 


This equation also gives the impedance Z, at the input end of a simple conduit 
of length (/;+ /2:), assuming Z,= 0 as before. From these considerations, it 
may be concluded that the filter system, with terminal conduits of length J; 
and /, should be as closely matched to free space as a simple conduit of length 
(1;+ /2), and that comparison with such a simple conduit ought to provide a 
useful measure of filter performance. It appears that the value of /; is not 
critical for such a comparison. 


A comparison of Equations (19) and (21) shows that in the first attenuation 
region, the curve for Za2 versus frequency can be made to coincide closely with 
the curve for Zo, for the case of the filter section of Fig. 1, (a), by a suitable 
choice of J. For the filter section of Fig. 1, (6), the curves for Za2 and Zp inter- 
sect near the center of the attenuation band. Hence, in this region, it ought to 
be possible to obtain impedance matching. The curves of Za2 for the terminal 
conduits used here are included in Fig. 3. 


To connect two filters of different characteristics in series, a coupling conduit 
may be required. If the filters have the same cross sections, they will have the 
same characteristic impedance over a large part of the first transmission band, 
and therefore they may well be connected directly. 


Experimental Method 


A method of measuring filter characteristics has already been suggested in 
the theory. The experimental arrangement is shown in Fig. 6. Readings of 


FILTER 
ouTPuT 
‘i 
L ] 


Fic. 6. Experimental arrangement. Sound source, Western Electric #555 receiver. 
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input and output sound levels were taken with the filter terminated in con- 
duits of length /, and /:, and with the filter replaced by a simple conduit of 
length (/;+ /2). The difference in attenuation in the two cases is the atten- 
uation due to the filter. 


Experimental Results 


Measurements have been carried out on several filters, of both circular and 
rectangular cross section. Two triple section filters, made up of cylindrical brass 


(b) 


L, cm. 


| 
| 


| 
a;,cm. | da, cm. 





30.0 | _— | 3.8 | —~0.685 
60.0 5; 5.0 —0.805 





60.0 1.65 +0.805 


30.0 | 1.65 | +0. 685 
i | 





Fic. 7. Triple section filter. (a)m<1. (b) m> 1. 


tubing, and designated #6a, #6b, and #7a, #7), are shown in Fig. 7. If b = — 8, 
then the attenuation characteristic of a three section filter is given by 


Pout = Pine* , (23) 


where, of course, 8 varies with frequency, Equation (16), and is a maximum 
at the center of the attenuation band. The terminal conduits are made of 
tubing 2.4 cm. in diameter, for which Z = poc/S = 9.3 acoustic ohms. Curves 
for the characteristic impedance Zo are shown in Fig. 3. Also shown are the 
output conduit impedances Z,2. The theoretical frequency limits of the atten- 
uation bands of each filter are indicated. The lengths of output conduit corres- 
ponding to the impedance curves of Fig. 3 are 3.3 cm. and 6.7 cm. for filters 
#6 and #7 respectively. Measurements were carried out using several lengths 
of output conduit, and the above values were found to be optimum. The 
attenuation curves obtained are shown in Figs. 8 and 9. The two filters were 
also connected in series, with #7 at the input end. The attenuation curve is 
shown in Fig. 10. 
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Fic. 8. Attenuation, filters #6a and #6b. 
1. Theoretical attenuation. 
2. Experimental attenuation, 1,= 3.3 cm., l2= 1.5 cm. 
3. Experimental attenuation, 1,;= 3.3 cm., la= 3.3 cm. 
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Fic. 9. Attenuation, filters #7a and #7b. 


1. Theoretical attenuation. 
2. Experimental attenuation, l,= 6.7 cm., la= 6.7 cm. 
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Fic. 10. Attenuation, filters #6a and #7a, #6b and #7b in series, #7 at input end. 
1. Theoretical attenuation, #7a, #7b. 
2. Theoretical attenuation, #6a, #6b. 
3. Experimental attenuation, 1, = 6.7 cm., le= 1.5 cm. 
4. Experimental attenuation, l,= 6.7 cm., l2= 3.3 cm. 


Attenuation measurements were carried out on a series of filters of rect- 
angular cross section, made up of both types of section shown in Fig. 1, and, as 
before, agreement with theory was satisfactory. In these filters, the expansions 
and contractions were confined to only one transverse dimension, the other 
remaining constant over the entire length of the filter. 


Conclusions 


It may be concluded from experimental measurements on filters designed 
according to the preceding theory that, up to frequencies to which measure- 
ments were carried, agreement with theory was satisfactory. 


The filters of rectangular cross section are particularly suitable for the con- 
struction of a filter unit of large cross section, since the two types of section 
can be placed side by side, with a single partition forming one boundary of both. 
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ON THE FLOW OF GASES AND WATER VAPOR THROUGH WOOD! 
By P. M. PFALZNER 


Abstract 


The quantities of water vapor passing through wood have been measured 
under the conditions of hydrodynamic flow and of kinetic diffusion. Air and» 
oxygen pass through wood under hydrodynamic flow conditions in accordance 
with the Knudsen-Poiseuille law. Water vapor passes through wood under 
hydrodynamic flow conditions in much larger quantities than those predicted 
by the law. The rate of passage of water vapor under diffusion conditions fol- 
lows the Fick diffusion law but with a diffusion coefficient that increases with 
relative humidity. The distinct features of the diffusion and hydrodynamic 
flow processes are pointed out. Porous, hygroscopic media permit the transfer 
of water vapor either by diffusion or by hydrodynamic flow, as determined by 
two distinct sets of conditions, while compact, hygroscopic media permit the ~ 
transfer of water vapor by diffusion only, regardless of the external conditions. 


1. Introduction 


In recent years it has become increasingly evident that the movement of 
water vapor through hygroscopic media is a complex phenomenon. Water 
vapor can move through such media more readily than can permanent gases 
such as oxygen, nitrogen, or even hydrogen. Any simple theory of the mech- 
anism by which the passage of gases and vapors takes place leads one to expect 
that small and light molecules will penetrate a porous body more easily than 
large and heavy molecules. The ease of penetration of porous hygroscopic media 


by water vapor points to some special mechanism facilitating the movement of 
water molecules. This mechanism must be sought in the fact that water vapor 
and other vapors are strongly sorbed by porous, hygroscopic bodies while the 
permanent gases are sorbed in negligible quantities only. 


Several mechanisms by which water molecules might be enabled to move 
more readily through porous, hygroscopic media have been suggested by a 
number of investigators. So far these have answered the problem in a quali- 
tative way only; the precise laws governing the movement of water vapor 
through porous bodies remain to be formulated. 


The simplest considerations suggest at least three mechanisms by a com- 
bination of which water vapor may move through a hygroscopic medium: 
(1) Since water vapor is a gas (close to its condensation point) it should pass 
through porous bodies in the same manner in which permanent gases do. The 
laws of kinetic diffusion of gases through porous interfaces have been well 
established; knowing the coefficient of diffusion for one permanent gas, one 
readily calculates the coefficients for other gases and, in particular, that for 
water vapor. Likewise, there exists a semi-empirical law of Knudsen—Poiseuille 
for the hydrodynamic or pressure flow of gases through porous bodies; again, 

1 Manuscript received January 20, 1950. 
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one may calculate the coefficients for the flow of water vapor from the coeffi- 
cients found for a permanent gas. In this way, one obtains the amount of water 
vapor that would be transported through a medium either by kinetic diffusion 
or by hydrodynamic flow, on the assumption that water vapor behaves—at 
least in part—as any other permanent gas. As pointed out already, this cal- 
culated amount is much less than the actually observed transfer of water 
vapor, and hence, to account for this excess transfer, there must exist additional 
mechanisms based on the hygroscopicity of the medium. 


(2) Water vapor can move through a porous body by capillary action. If the 
capillaries become filled, the water will tend to rise in accordance with well 
known laws. Water is certainly drawn into porous bodies in this way, but in 
general such movement is confined to conditions in which the moisture content 
of the medium is quite close to the saturation value. It is known from experi- 
mental evidence, however, that an excess transfer of water vapor takes place 
even at the very lowest moisture contents, and hence capillary action cannot 
wholly explain the observed facts. 

(3) Evidence from a variety of experiments has indicated that adsorbed 
water molecules are not static but are free to move over the surface of the 
adsorbing medium. This suggests that water vapor may move through the 
medium as an adsorbed layer on the internal surfaces. It is this mechanism 
that presents the greatest interest to research. 


The experiments to be described below are a first attempt at separating the 
contributions of mechanisms (1) and (3) to the vapor flow. In these experi- 
ments, the transfer of gases has been studied under hydrodynamic flow con- 
ditions rather than diffusion. Diffusion experiments with permanent gases 
present considerable difficulties while hydrodynamic flow measurements are 
relatively simple. 


It should be pointed out here that there exists a good deal of confusion in 
the literature in the use of the term ‘‘diffusion’’ (which is often applied where 
the term “hydrodynamic flow” is more appropriate). The term “hydrodynamic 
flow”’ is here used to describe fluid flow under a pressure head, while the term 
“‘diffusion’’ of gases should be restricted to describe the exchange of gases from 
opposite sides of a barrier (real or imaginary) or a porous interface, where there 
is a partial pressure gradient but where the total pressure on both sides remains 
constant throughout. Such exchange of gases will occur, for instance, in the walls 
of a building where the total pressure is atmospheric both inside and outside but 
where there is a difference in water vapor pressure between the exterior and 
interior of the walls. The water vapor will accordingly diffuse through the wall 
until equilibrium is established. At first sight, it would appear that such a 
process should obey the same laws as hydrodynamic flow, since it is known 
from kinetic gas theory that within certain limits in a mixture of gases each 
gas behaves as if it alone were present. Consequently, one might expect that 
the presence of air on both sides of an interface ought not to influence the 
movement of water vapor across this interface. This is not so, however. Not 
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only does there also take place a counter movement of an air mass equal to 
the mass of the diffusing vapor, but the vapor also has to diffuse through the 
air towards and away from the interface. The counter movement of air will 
be very small since the mass of water vapor in the air at normal temperatures 
is quite small—at 20°C., 100% relative humidity, the absolute humidity is 
17.5 gm. of water vapor per cubic meter of air, or the mole fraction of water 
vapor is 0.0218. Such a counter movement of the air, however, does not take 
place when the vapor diffuses from one side of an interface to the other and 
condenses there. No air will then flow in a direction opposite to that of the 
vapor flow since the condensing vapor does not increase the total pressure. 















The fundamental law of diffusion of two gases is, according to Stefan (16), 
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Pe = —— 
CP dx 
where p = density of the gas 
u = velocity of the gas 
RT ; 
k= 7, = a constant for any given temperature and gas 
P = total pressure 
= pit pr 
p = partial pressure 






C = coefficient of resistance. 





The subscripts refer to the two diffusing gases. 






If only one gas is diffusing, the equation governing this movement is 
aks don 
C(P — p:) dx 






(3) 


ay = = 






The hydrodynamic flow of a gas through a circular tube is given by Poise- 


uille’s law i le 
5 2 ’ (4) 








Gp = K(fi — po) 






where G = volume of gas flowing in unit time, 
. ert airs ; } . 
K + = the Poiseuille constant for a circular tube, of radius 7, 
8nc 





length d, and gas of viscosity 7 
p, and p2= the pressures at the two ends of the tube, and 






p = pressure at which the volume G is measured. 






In the experiments described below, the procedure was: 
(a) To measure the flow of a permanent gas through a porous body, 
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(b) To calculate the flow to be expected for water vapor under the same 
conditions, 


(c) To measure the flow of water vapor and compare with the calculated 
value as obtained from (0). 


As a simple concept, it is then assumed that any excess of observed over 
calculated amount of water vapor flow is due to either or both of the mech- 
anisms (2) and (3) described above. 


In addition to the hydrodynamic flow measurements with both water vapor 
and permanent gases (oxygen and air were used), diffusion experiments were 
also made with water vapor only. These will illustrate the difference between 
hydrodynamic flow and diffusion of water vapor. 


The hygroscopic medium chosen for these experiments was wood. This was 
done for a number of practical reasons. Wood has a fairly uniform pore struc- 
ture with capillaries ranging from coarse to very fine (14). It adsorbs water 
vapor strongly but does not sorb air or oxygen to any extent. Much work has 
already been done by various investigators on the diffusion of water vapor 
through wood and charcoal but not on the hydrodynamic flow of water vapor. 


2. Experiments with Hydrodynamic Flow of Air, Oxygen, and 
Water Vapor 


(a) Apparatus 


The apparatus was designed to permit observations to be made at constant 
pressures so that a stationary concentration gradient would be maintained 
throughout the disk of wood. As shown in Fig. 1, there are two absolute mer- 
cury manometers M, and M, with attached movable mercury reservoirs R; 
and R, on each side of the cell A in which the disk is clamped. In the figure, 
gas flow takes place from Ms, the high pressure or inflow side, towards M, the 
low pressure or outflow side, either upwards or downwards through the disk, 
as determined by the position of the threeway stopcocks G, G. The mercury 
level in M, is held just below the tungsten contact F; whenever the pressure 
falls owing to the flow of gas across the disk, contact is made at F and a relay 
closes a circuit which causes gas to be generated electrolytically in the bulb 
Bz filled with water containing a little sulphuric acid. The pressure increase in 
Bz causes the mercury column in the burette D,» to rise until an amount equal 
to the volume of gas flowing through the disk has been displaced, and the circuit 
is once more opened at F. The rise of the mercury column in D,, measured with 
a cathetometer, serves as a measure of the rate of flow of gas across the disk. 
At the low pressure side a similar relay is actuated by the tungsten contact S 
sealed into manometer /, and starts a pump acting on the bulb B,; which 
causes the mercury in D, to fall. As soon as the contact S opens, the pump 
stops and air is permitted to leak back through the tap H into the bulb By. 
This ensures continual opening and closing of the circuit. The pull of the 
pump may be suitably throttled by adjusting taps H and J. At the lowest 
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Fic. 1. Schematic drawing of apparatus for measuring the flow of gases through wood. 


pressures, when the flow takes place into vacuum, the pump can be switched 
to act directly on the outflow side of the disk. 


A differential manometer M was added during the course of the experiments 
to obtain direct readings of the pressure difference maintained across the disk. 
It was then unnecessary to read manometer M,, as the absolute pressure was 
read on M, and the pressure head on M; the latter carried the tungsten con- 
tacts controlling the manostat pump. The same pressure head could be obtained 
accurately at all mean pressures as long as the amount of mercury in M was 
kept unchanged. 


The disk of wood was sealed at its edge with several coats of varnish and 
beeswax, the latter applied while the disk was clamped in position with the 
apparatus under vacuum. The metal flanges holding the disk were faced with 
rubber washers and sealed to the glass tubing with de Khotinsky wax. 


(b) Method of Taking Readings 


‘The apparatus was exhausted with a mercury diffusion pump and tested for 
leaks. Gas was admitted through Z until the desired pressure was obtained at 
M;,. The stopcocks G were then adjusted for upward or downward flow, and 
additional gas was admitted until the desired pressure head was reached. The 
mercury levels in M,; and M:2 were adjusted so as to be just below the contacts 
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Sand F respectively. (After the manometer M had been installed, the pressure 
was increased by adding more gas until the mercury in M was just below the 
contact K.) No measurements were taken during the first minutes of flow, in 
order to permit establishment of a stationary pressure gradient inside the wood. 
The time was then measured for a volume of from 30 to 50 cc. of gas to pass 
through the disk as indicated by the rise of the mercury column in the cali- 
brated burette D2, the volume readings being taken with a cathetometer, the 
time with a stopwatch. The room temperature, which didnot vary more than 
1°C. during any one setting of the pressure head and mean pressure, and not 
more than 2°C. over a complete set of measurements for one gas, was read at 
the beginning and end of each observation. (See Appendix II for effect of 
temperature changes.) For each mean pressure two readings were taken, the 
first with the flow upwards through the disk, the second downwards. 


Both the air and the oxygen were dried over magnesium perchlorate before 
being admitted into the apparatus. The oxygen was taken from a pressure 
tank of the Dominion Oxygen Co. 


The same circular disk of wood was used throughout the experiment. It was 
a piece of spruce wood cut from a radial section, i.e., the flow took place in a 
direction tangential to the yearly rings. Size and other characteristics as deter- 
mined by methods described in Appendix I were as follows: 


Thickness of disk at 0% humidity 0.629 cm. 


Diameter of disk 3.39 cm. 

Area of disk available for flow 9.02 sq. cm. 
Void fraction when dry, Wo 0.78 cc. per cc. 
Dry bulk density, pp 0.33 gm. per cc. 
Density of wood substance, p 1.48 gm. per cc. 


(c) Theory 
(1) The Streamline Flow of Gases Through a Single Capillary 
Poiseuille’s law is, as stated above, 
Gp = K(pi— pa) EP, (4) 


with the symbols defined as before. 


(2) The Flow of Rarefied Gases 


At pressures where the mean free path of the gas molecules is of the same 
order of magnitude, or larger, than the diameter of the capillary, the flow no 
longer obeys Poiseuille’s law. At these pressures the flow becomes much larger 
than that law indicates. Knudsen (10), on the basis of kinetic theory and with 
the added assumption that a gas molecule may be reflected in any direction 
from the wall of a capillary regardless of the angle of incidence, developed an 
equation of the form 


Gp = b(bi— p2), (5) 
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where D 


gas constant, 

absolute temperature, 

molecular weight of gas, 

radius of capillary, 

length of capillary. 

This equation has been shown to be valid at very low pressures (7), the 

volume flowing being inversely proportional to the square root of the molecular 
weight and directly proportional to the third power of the radius. The type 


of flow occurring under these low pressures has been termed ‘‘molecular flow” 
to distinguish it from Poiseuille or ‘‘streamline flow’’. 


In the pressure region intermediate to molecular and Poiseuille flow, Knudsen 
combined Equations (4) and (5) with a variable factor y obtaining 


i -bi t+ peo 
Gp = K~—~ (pi — bs) + y(bi — Pa). (6) 


Knudsen found that y = i GP 
1+ C.P 


where C, = 2.00 res , Ce = 2.47 4/ BE s 


(7) 


M M 2 


P = mean pressure. 


Equation (6) is valid for a wide pressure range; at low pressures y — 1, and 
the first term in (6) becomes negligible, thus giving Knudsen’s law of molecular 
; Ci ; 
flow, while at high pressures y —> —, and the equation reduces to the Poise- 
2 ( f 
° . « iit . 1p 
uille law of streamline flow with an added constant. The curve of 


Pi— P 


°c ° 1 yy) ana : 
specific flow, against mean pressure —— -_ has a minimum ata pressure \ 
2 


the mean free path of the gas molecules is approximately five times the « 
lary radius. 


It must be emphasized that Equation (6) takes no account of sorption, and 
in fact Knudsen’s derivation assumes a uniform concentration of the gas mole- 
cules in any cross-section of the capillary. 


(3) The Flow of Gases Through Porous Systems 


‘When a gas flows through capillaries whose radii are of the order of magni- 
tude of the mean free path of the gas molecules at normal pressures, the con- 
ditions are formally the same as for rarefied gases flowing through coarse 
capillaries as described under (2) above. While Knudsen developed his Equa- 
tion (6) for rarefied gases, Adzumi (2) has shown that an equation of the same 
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form, with y replaced by a constant, is applicable to gases flowing through a 
porous system. Any complex system of capillaries will be made up of capillaries 
of a wide range of radii, some so narrow as to fall in the region where molecular 
flow occurs at normal pressures, and others wide enough to permit streamline 
flow. On the basis of his experiments, Adzumi replaced y by a constant = 0.9, 
and wrote the equation 


Gp = (ZeP +09 x 4/3 x Vie 4/ ETP) (ps — ps (9) 
n 4 


where E and F are functions of the radii and lengths of the capillaries. The 
form of these functions is, of course, unknown owing to the complex structure 
of most capillary systems. The values of E and F, however, may be obtained 
from a graph of specific flow Gp/(p1— p2) against mean pressure P which will 
be a straight line (if the equation correctly describes the flow) having a slope 


of rE/8n and an intercept of 0.9 XK 4/3 X V/ 25 - F. 


(4) Comparison of Flow of Different Gases Through the Same Porous System 


If two gases of molecular weights M, and Mz, and viscosities 7; and 72 flow 
through the same capillary system under the same conditions, then from (9) 


Gip = (AiP + B,)(pi— po), (10) 
and = 

Gop = (A2P + Bz)(pi-— po), (11) 
where A us sk (12) 

8n a 

B = 09X4/3 X~W/2n eS F, (13) 

and hence A, om (14) 
2 71 

B,_ , /M, 

and B, = ft: (15) 


Thus if A and B are found for any one gas, the corresponding constants for any 
other gas may be computed from Equations (14) and (15). 
(d) Results for Hydrodynamic Flow of Air and Oxygen 
(1) Air 
The volume of dry air flowing at mean pressures from 4.4 cm. to 80.0 cm. of 
mercury with a constant pressure head of 5.50 cm. of mercury was measured. 
In Fig. 2 values of the specific flow Gp/(pi— 2) are plotted against the mean 
pressure (pi+ p2)/2. The slopes and intercepts taken from the graph are given 
in Table I. 
(2) Oxygen 
The flow of oxygen was similarly measured, and the slope and intercept 
obtained from the graph for oxygen (Fig. 2) are also given in Table I, together 
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Fic. 2. Flow rates of air and oxygen through wood at 20.8°C. 










with the values calculated from the measured air flow by use of Equations 
(14) and (15). The calculated and observed figures for oxygen flow based on 
air-flow data are within 3% of each other. Other investigators, such as Tom- 
linson and Flood (17), have recently confirmed that the Knudsen—Poiseuille 
equation accurately predicts the flow rates of permanent gases. They used a 
method similar to the one described here, except that their medium of flow 
was activated carbon. Taking the flow of nitrogen as basis, they found good 
agreement between calculated and observed flow rates of carbon dioxide and 
helium. For hydrogen they found a somewhat lower observed rate than that 
calculated from the equation. 












The values for the viscosities and molecular weights of air, oxygen, and water 
vapor at 20.8°C. which were used in the calculations were calculated from 
values given in Landolt—Bérnstein Tabellen (10). These figures are given in 
Table II, together with the viscosity of air containing water vapor corres- 







TABLE I 
HYDRODYNAMIC FLOW OF AIR AND OXYGEN 

















A (slope) B (intercept) 








Cc./sec./mm. Hg Cc./sec. 


Se ethos Oe ad 0.000139 0.0400 
MNO ok occ are eens 0.000126 0.0393 
Oxygen (calc.).... .| 0.000123 0.0380 


% Difference...... 3% 
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ponding to 75% relative humidity at 20.8°C., calculated from a relation for the 
viscosity of gas mixtures given by Adzumi (1). 


TABLE II 


VISCOSITY AND MOLECULAR WEIGHT 
OF AIR, OXYGEN, AND WATER VAPOR 


Viscosity X 10* Molecular weight 





c.g.s. units 
Far tary). on «2 32s o 1.814 29 
Air (75% R.H.)..... 1.806 29 
CURE ce isdcccves 2.054 32 
Water vapor........ 0.949 18 


(3) Air at 75% R.H. 

The flow of air containing water vapor corresponding to 75% relative 
humidity at 20.8°C. was measured in order to discover whether the presence 
of water vapor in the air had any influence on the rate at which the air passed 
through the wood. The presence of water vapor has only a slight influence on 
the viscosity of the mixture. As seen from Table II, for air.at 20.8°C., 75% 
R.H. the change is 0.45%. Hence, there should be no noticeable change in the 
measured air volume if water vapor behaves just like any permanent gas. 
However, owing to adsorption and swelling of the wood, changes in the flow 
rate might be expected. 


In order to obtain 75% relative humidity, air was bubbled through saturated 
sodium chloride solution and circulated through the apparatus by means of a 
reciprocating pump, past wet and dry bulb thermometers sealed into the 
system, and back to the saturated salt solution. After the thermometers had 
been registering a steady relative humidity of 75% for several hours, the circu- 
lation was halted, and the flow through the wood was measured in the manner 
described above for air and oxygen. It was found that the flow was not meas- 
urably different from that at 0% humidity. Measurements with the same 
negative result were obtained after the wood had been exposed to the humid 
atmosphere for more than a week. 


The reason—apart from the negligible change in viscosity as pointed out 
above—is that the volume of moisture contained in the air and in the wood 
under the conditions of the experiment is small. The absolute humidity at 
75% R.H., 20.8°C. is 12.5 gm. per cu. m., or 0.0156 mole fraction H,O in air; 
the moisture content of the wood at 20°C. is approximately 14% of the dry 
weight (9, p. 57), and for wood of a bulk density of 0.3 gm. per cc. there will 
be approximately 0.04 gm. of water content per cubic centimeter. Hence, if 
the density of the sorbed water is not less than that of free water (the opposite, 
higher density of the sorbed water than that of free water, has been found by 
various investigators, see Stamm, 15, p. 58), about 4% or less of the volume 
of the wood will be occupied by water, and the void fraction will be reduced 
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by this amount. A 4% reduction in the void fraction will not appreciably 
reduce the volume available for flow, especially since any condensation of the 
water will occur initially in the narrowest capillaries which do not contribute 
measurably to the flow. 







(e) The Hydrodynamic Flow of Water Vapor 


(1) Method 

In order to measure the flow of water vapor through the wood a gravimetric 
method was adopted for the following reasons. The flow takes place at very 
low pressures (the saturation pressure of water vapor at 20°C. is 17.54 mm. of 
mercury) and the volume flowing is therefore very small. Hence, it is more 
convenient to measure the weight of water vapor passing through the disk. 
To obtain a given vapor pressure, a flask containing a suitable saturated salt 
solution (see Table IV) was sealed to the apparatus and kept in a water bath 
thermostatically controlled at 20°C. Distilled water was used to make up the 
saturated salt solutions, and distilled water alone was used to give the highest 
relative humidity, 100%, at 20°C. It must be realized, however, that the rela- 
tive humidity at the disk of wood was at all times slightly lower than that 
given by the saturated salt solution at 20°C. since the temperature of the 
apparatus and wood was always slightly above 20°C. to avoid condensation 
of the water vapor. To ensure saturation of the salt solutions at all times, they 
were agitated by means of a magnetic stirrer. The air was exhausted from the 
apparatus. The flow then took place through the disk into vacuum maintained 
by continuous pumping. Between the pump and a Schwartz drying tube con- 
taining magnesium perchlorate which served to collect the outflowing water 
vapor, a calcium chloride drying tower was inserted to prevent moisture from 
the pump rising into the Schwartz tube. The vapor pressure was read directly 
on the mercury manometer M (Fig. 1) with a cathetometer. The Schwartz 
tube was taken down and weighed at intervals of about 20 hrs. The readings 
were repeated at the same vapor pressure at least four times, or longer, if 
necessary, until the change in weights recorded remained constant. 

































(2) Results 


From the increase in weight of the Schwartz tube, the following values 
(Table III) were calculated: amount of water vapor passing through the 
disk in 


(a) grams per hour, Wo, 










(b) grams per hour per millimeter mercury pressure head, W,, 


(c) specific flow, Gp/(pi— pe). 












In Fig. 3, the specific flow of water vapor as measured is plotted against the 
mean pressure. In addition, the curve predicted by the Knudsen—Poiseuille 
equation for gas flow has been calculated for water vapor from the air-flow 
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data and drawn on the same plot. The slope A and intercept B of this curve 
were calculated by use of Equations (14) and (15) and are 


A = 0.000265 cc. per sec. per mm. mercury 
B = 0.0507 cc. per- sec. 


TABLE III 
HYDRODYNAMIC FLOW OF WATER VAPOR 








bi + po (a) ) (c) 
i- 2) | a 4 
aS fis 2 Wo gm./hr. W, gm./hr./mm. Hg ma 9 8 b:) 
mm. Hg xX 10 X 104 x 10? 
2.20 1.10 10.60 4.82 13.60 
5.58 2.79 16.08 2.88 8.13 
7.14 3.57 19.64 2.75 7.76 
8.13 4.06 22.35 2.75 7.76 
9.43 4.72 23.12 2.45 6.91 
11.58 5.79 25.63 2.21 6.23 
13.83 6.92 29.42 2.13 6.01 
14.91 7.46 32.37 2.17 6.12 
16.00 8.00 50.75 3.17 8.94 
17.83 8.92 47.31 2.65 7.48 
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Fic. 3. Flow rates of water vapor through wood at 20°C. 







It will be seen that the experimental curve is totally different from the 
calculated curve. At the lowest mean pressures, the observed value of the 
specific flow is almost three times that given by the Knudsen—Poiseuille equa- 
tion. The specific flow decreases rapidly to a minimum at a mean pressure of 
about 6.5 mm. of mercury, and then rises once more as saturation pressure is 
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approached. The Knudsen—Poiseuille curve, which lies below the experimental 
values is, of course, a straight line. It must be noted that the experimental 
curve was obtained from the average flow rates over the range extending from 
zero pressure to that given by the saturated salt solutions, corrected for unit 
pressure difference (integrated flow rate). The Knudsen—Poiseuille curves for 
air and oxygen were obtained as average flow rates over a constant range of 
pressure difference for varying mean pressures, and were also corrected for 
unit pressure difference (differential flow rate). Since the Knudsen—Poiseuille 
curves are straight lines, a straight line should always be obtained for gases or 
vapors obeying the Knudsen—Poiseuille flow relation regardless of whether the 
integrated or differential flow rates are plotted. The differential flow rates for 
water vapor would yield a curve of a shape different from that of the integrated 
curve plotted in Fig. 3, but it would certainly not be a straight line. It probably 
has a more pronounced minimum than the curve drawn in Fig. 3. 

















Tomlinson and Flood (17) have measured ‘the flow rates for the vapors 
diethyl ether and ethyl chloride through activated carbon rods. They found 
that the flow rates of these vapors, which are adsorbed, are considerably larger 
than the Knudsen-Poiseuille relation predicts, and their flow rates for diethyl 
ether show a minimum at low pressures similar to that obtained here for water 
vapor. The suggestion has been made by Tomlinson and Flood that the mole- 
cules in the adsorbed layer are mobile and that the excess flow is determined 
by the concentration gradient in the adsorbed layer. Owing to the S-shape 
of the adsorption isotherm of water vapor on wood, the concentration gradient 
for equal pressure differences will be large at very low pressures and near 
saturation and, hence, the flow will be large in these regions. The curve of 
Fig. 3, which is that of the integrated flow reduced to unit pressure difference, 
supports this view. If the differential curve had been obtained, the minimum 
would be more pronounced, as pointed out above. 

















According to Pollard and Present (12), the expression for the interdiffusion 
of two gases in long capillary tubes reduces to the Knudsen equation for mean 
free paths much greater than the capillary radius; according to these same 
investigators, the curve of specific flow against mean pressure will show a 
minimum at low pressures only if the flow takes place through uniform capil- 
laries; in porous media, owing to the effect produced by the range of capillary 
radii, there will be no minimum. This seems to hold true for nonadsorbed gases 
such as oxygen and air investigated here, and hydrogen, helium, carbon dioxide, 
and nitrogen as measured by Tomlinson and Flood. It does not hold for water 
vapor as shown here, or diethyl ether as pointed out above. 















3. Experiments on Diffusion of Water Vapour 


(a) Method 


In order to measure the rate of diffusion of water vapor through wood, a 
dry-cell method similar to that employed by Babbitt (3) in measuring the per- 
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meability of building papers was adopted. Circular disks of spruce wood of 
three different thicknesses (0.376, 0.624, and 1.276 cm.) with an area of 37 
sq. cm. were sealed across the mouth of aluminum dishes containing calcium 
chloride. The wood was taken from the same board from which the disk used 
in the hydrodynamic flow measurements had been cut. The aluminum dishes 
were exposed to humid air in a constant humidity chamber kept at 91°F. 
(32.8°C.). This temperature, rather than 20°C. as in the hydrodynamic flow 
experiments, was chosen here because it was simpler to control thermostati- 
cally a closed chamber at a temperature above room temperature than below 
it. Salt solutions were used as before to establish constant relative humidities. 
These are given in Table IV for both 20°C. and 32.8°C. 


TABLE IV 


SATURATED SALT SOLUTIONS TO ESTABLISH 
DEFINITE WATER VAPOR PRESSURES 














————— a 
20°C.* 32.8°C.f 

Salt Vapor pressure, RA: Vapor pressure, R.H. 
mm. Hg. % mm. Hg. % 
Nass 6. ah. on! ciare 16.0 91 — - 
Orcs cron swe 14.9 85 32.8 88 
SMS viata choos 5 13.8 79 28.0 75 
PNM IN OR i ceieyecniets | 11.6 66 —_ = 
Ca(NOs3)2 ete € mae 9.4 54 aie = 
LOC be SN eae 8.1 46 _ - 
MSE Be os isia ie toe 41 12.3 32 
CH;COOK...... 5.6 32 — - 
Mas iis ca cee oe 2.2 12.5 — - 











* The vapor pressures at 20°C. were measured directly on manometer M for each flow experiment 
and averaged, and the corresponding relative humidities calculated. These are not in every case 
identical with the ones given in Landolt-Boérnstein Tables. 

+ The vapor pressures at 32.8°C. were obtained by plotting data given by O’Brien (11). The cor- 
responding relative humidities were calculated. 


(b) Theory 

The diffusion equation of Stefan has been given above. For only one gas 
diffusing, it is 
_ hits dn, @) 
C(P — pi) dx 
where the symbols have been defined previously. This equation can be inte- 
grated as follows: Substituting kip,1= pi, we obtain 
U1p\dx = ke se 
C (P — pi) 


which gives, on integration with respect to x and p, 


Aawi=— 


upid = = log (P — pi) +A, 


where d = thickness of the medium, 
and A = aconstant of integration. 
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Now for ~1 = 0, the condition is that u,= 0, hence 
ko 
A = ——log P 
C ig 


Ppt. 
P 





: ke 
. tind = —lo 
1P1 C g 


velocity times density, 


now U1 p1 
mass flow per unit time per unit area of cross section, 
W 
At 


where W = mass flowing in time ¢ across an area A; 
i.e., the equation can be written 


Wd Bnd ko P = Pi 


— Sige (16) 











= —log ie. 
At C P 
The common form of Fick’s diffusion law is 
Wd 
— = DAp, 
At ? 
where Ap = p1— pz», the partial pressure gradient. (17) 
If p2= 0 at all times, Ap = fi, 
; Wd 
1.€., = Dp. 18 
7 pr (18) 


Equation (16) can be expanded into 


s--25+36)-26)+36)—-1 
At - BBS P 3 \P "7 P oy 


which is, on a first approximation, 





We wee, (19) 


Al CF 


Hence D, the diffusion coefficient in (18), is equal to — a approximately, 


while the value of the true diffusion coefficient is given by 


Wd _—, P-p, 
= Dloe—_"!, (16) 
‘>; j 


where D= * (20) 
Cc 
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(c) Results 


The ratio= in Equation (18) was calculated from the steady rate of increase 
C 


in weight of the aluminum dishes. It is given in Table V in grams per 24 hr. per 
square meter per mm. of mercury and is plotted against mean relative humidity 
in Fig. 4. It is seen that the permeability increases rapidly with mean humidity 
and that it depends on the thickness of the medium. The permeability decreases 
steadily at decreasing humidity, as has also been shown by Babbitt (4) and 
King (6). The curve is again that of the integrated water vapor transfer, i.e., 
the mean rate of diffusion over the range of pressures extending from zero to 


Hg' 
NV 
°o 





Te 0.376 Cm. _ 
or 0.624Cm. 

N 1.276 Cm, 

E hi 

ee \ 
Alo 


° 5 10 15 20 25 30 35 40 45 50 
MEAN RELATIVE HUMIDITY,% 


Fic. 4. Diffusion rates of water vapor through wood at 32.8°C. 


that maintained in the constant humidity chamber. The coefficients are 
plotted against the mean relative humidity rather than the mean vapor pres- 
sure, since the moisture content is a function of the relative rather than 
absolute humidity and since media of the same moisture content rather than 
at the same mean vapor pressure are comparable with regard to permeability. 
The specific permeability or coefficient of diffusion, D, is also given in Table V, 
































TABLE V 
DIFFUSION OF WATER VAPOR 
Thickness, cm. 
0.376 | 0.624 | 1.276 | 0.376 | 0.624 | 1.276 
Pp, Mean - 
mm. Hg. | R.H., | Wa D, 
% | D/d = W/tAp,, a av. 
| 
gm./24 hr./sq. m./mm. Hg. | gm. cm./24 hr./sq. m./mm.Hg. 
ne { 
12.3 16.5 | 0.53 0.31 0.13 0.20 0.20 0.17 0.19 
28.0 37.5 1.39 0.80 0.55 0.52 0.50 0.70 0.57 
32.8 44.0 1.98 0.98 0.64 0.74 0.61 0.81 0.72 
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and it will be noted that for the same humidity this value is independent of 
the thickness of the medium, within the experimental errors. 







c The true coefficient of diffusion, D, has been calculated in Table VI and is 
‘ plotted against mean relative humidity in Fig. 5. The ratio D/D is given in 
y 9r- 

y 

5 
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Fic. 5. Diffusion rates of water vapor through wood. Du from hydrodynamic flow; D from 
diffusion. 


O 
wo 












the table; it is seen to be nearly constant and has an average value of 75.6; 
this should follow also from the fact that D = k2/C while D = k2/CP approxi- 
mately, and P = 76.0 cm. of mercury. The approximation made in Equation 
(19) is seen to be justified, i.e., Equation (18) may be used to calculate relative 
values of the diffusion coefficients. 








It is not possible to obtain a value corresponding to D from hydrodynamic 
flow experiments since there is no equivalent for P, the total pressure. It is 










TABLE VI 
THE DIFFUSION COEFFICIENT D FOR WATER VAPOR 























mean R. H. P—-p D= ml r—fi D 
ian ‘Hg. = mm. Hg. ni" FP D/D 
| 16.5 747.7 14.5 76.3 
37.5 732.0 42.6 74.7 
44.0 727.8 54.5 75.7 
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possible, however, to calculate a value which is here called Dy corresponding 
to D from these data. This has been done in Table VII where Dy has been 
calculated from the W, values of Table III. In Fig. 5 the average value of D from 
diffusion experiments (Table V) has been plotted against mean relative humid- 
ity together with the values of Dy from hydrodynamic flow experiments (Table 
VII). The two curves, both of which represent integrated vapor transfer, are 
entirely different. Quantitatively, the Dy values are about 10 times as large 
as the D values. D gradually approaches zero at low humidities while Dy passes 
through a minimum and then rises sharply. Thus besides the quantitative 
differences there are outstanding qualitative ones. The absence of air probably 
facilitates the passage of vapor through the medium and this may partly 
account for the quantitative difference, but there must be additional mech- 
anisms at work to account for the qualitative phenomenon of increased transfer 
of vapor under hydrodynamic flow conditions at low humidities as compared 
to the gradual approach to zero transfer under diffusion conditions. 


TABLE VII 


THE DIFFUSION COEFFICIENT DH FOR WATER VAPOR 
FROM HYDRODYNAMIC FLOW MEASUREMENTS 





x Corresponding Du= Wid X 24 X 104 
ace fe mean R.H., A 
F : 7% gm. cm./24 hr./sq. m./mm. Hg 
2.20 6.8 8.08 
5.58 16 4.83 
7.14 20.5 4.62 
8.13 23 4.62 
9.43 27 4.10 
11.58 33 3.70 
13.83 39.5 3.57 
14.91 42.5 3.63 
16.00 45.5 5.31 
17.83 50 4.44 








King (6) has measured the diffusion of water vapor through horn keratin and 
obtains a curve showing the variation of the integrated mass transfer with 
relative humidity. His curves and those of Fig. 4 (which represent the inte- 
grated mass transfer of water vapor through the wooden disks of varying thick- 
nesses plotted against mean relative humidity) have the same shape when 
reduced to the same conditions. King’s original curves represent values which 
have not been reduced to unit vapor pressure difference, and his abscissae are 
the relative humidities on the high pressure side in place of the mean relative 
humidity of the medium. Recalculating King’s figures, it is found that his 
values of D range from 0.03 at about 30% mean R.H. to 0.45 gm. cm. per 24 hr. 
per sq. m. per mm. of mercury at about 45% mean R.H. This compares with 
our values from about 0.40 at 30% mean R.H. to 0.72 gm. cm. per 24 hr. per 
sq. m. per mm. of mercury at about 44% mean R.H.; it is seen that our values 
are higher throughout but do not rise as sharply as King’s at higher humidities. 
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It follows that the diffusion coefficient for keratin changes more rapidly than 
that for wood as the humidity approaches saturation. 


It is interesting to note the conditions of King’s experiment. He excluded all 
air from both sides of a thin (5 X 107% cm.) keratin membrane and kept one 
side at zero pressure while the other was exposed to water vapor. These are 
the conditions we have so far termed those of hydrodynamic flow; nevertheless 
King applies the ordinary diffusion equation and obtains curves which as 
stated above are identical with our curves of Fig. 4 where the mass transfer 
due to what we have so far termed kinetic diffusion through a porous medium 
is plotted. In Fig. 5 we plotted data from hydrodynamic flow measurements 
using the diffusion equation to calculate Dy and we noted that the upper curve 
of Fig. 5 is totally different from the true diffusion curve. The reasons are to 
be sought in the structure of the flow medium. Keratin is not porous in the 
same sense as wood. Keratin may have a void fraction of the same order of 
total magnitude as some woods, but the voids are those of intermolecular spaces 
and not parts of a connected system of pores. Permanent gases will move 
through keratin—which may be termed a compact medium—at very low rates, 
if at all. Water vapor, however, is sorbed by keratin, swelling takes place, and 
the characteristic hysteresis loop of the absorption isotherm is observed exactly 
as for a porous medium such as wood. Water vapor will accordingly move in 
the adsorbed film by a kinetic process—but there will of course be no flow 
through capillaries. 


We can now elaborate the distinction previously drawn between diffusion 
and hydrodynamic flow. Diffusion of water vapor or any other adsorbable 
vapor takes place according to Equation (16), through either porous or com- 
pact media when there is an equal total pressure at the two faces of the medium 
with partial pressure gradients across the medium; or when there is a hydro- 
dynamic pressure gradient of a single vapor across a compact medium which 
sorbs the vapor. Hydrodynamic flow of a gas or a vapor takes place only 
through a porous medium, when there is a hydrodynamic pressure gradient 
of a single gas or vapor across the medium. The flow of nonadsorbable gases 
obeys Equation (9) while the flow of adsorbable vapors takes place in excess 
of the amount given by this equation. 


In other words, vapor transfer through a compact medium will take place 
under the conditions of true diffusion (i.e., in the presence of air with equal 
total pressure on both sides) as well as under the hydrodynamic conditions 
which King actually employed; in both cases the diffusion equation (16) is 
applicable, and similar curves will be obtained for the variation of the diffusion 
coefficient with relative humidity. King found that the presence of small 
amounts of air reduced the mass transfer of water vapor but he was able to 
account for this by the reduction in vapor pressure at the inflowing face of 
the keratin membrane which is caused by the resistance of the air through 
which the vapor had to diffuse when passing from the source to the medium. 
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King has successfully applied an analysis of the diffusion process given by 
Cassie (5) for wool to his experiments with keratin. We hope to discuss at a 
later date the application of the same analysis to a porous medium such as wood. 
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APPENDIX 


I. Preliminary Experiments 


The void fraction, density of wood substance, bulk density, and percentage 
swelling of the wood were determined. Small blocks of wood were dried to 
constant weight G in a vacuum desiccator at 104°C., then saturated under 
vacuum with water to weight G,, weighed while immersed in water, Gyy, and 
while immersed in benzene, Gy. The density p is then computed from 


_ ee f 
G = Guu 








where p, is the density of water at the temperature of the experiment. The 
void fraction W of the swollen wood is given by 
Gun + Vey > G 
Vey 
G 
Vp 








W = 


=1-— 





’ 
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where V is the volume of the swollen wood and may be measured or calculated 
from 


os Gun — Gus : 
Ppa PH 
where pz, is the density of benzene. The bulk density p, of the swollen wood is 
given by 
Pw = (1 — Wp. 


The figures for percentage swelling of the wood when taking up water from 
the oven-dry to the water-logged state were 0.8%, 1.5%, and 8.3% respectively 
for longitudinal, radial, and tangential swelling; the volumetric change was 
therefore about 10%. 


The bulk density pp of the dry wood is therefore given by 


p _- Fe 
”  -e 


and the void fraction Wo of the dry wood by 


The calculations gave the following results: 
W = 0.80 cc. per cc. 
Wo = 0.78 cc. per cc. 
1.48 gm. per cc. 
0.30 gm. per cc. 
0.33 gm. per cc. 


II. The Effect of Small Temperature Changes on the Constants 
of the Knudsen-Poiseuille Equation 


Equation (9) may be written 


fw AR Et. 
pPi-—pe 1 


; ‘ A ‘ 5. She 
The fractional change in slope — , for a change 6¢ in temperature, is given by 


and similarly, the fractional change in intercept B./T, for a change 6t in tem- 
perature, is 
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Now, the temperature dependence of the viscosity is given by the following 
laws: 

(a) For air: 1 = mo— a(23°— #) (10) 

(6) For oxygen: 7 = no(1 + dt)* (10) 

(c) For water vapor: » = no+ a(t — 100) (13). 











For air, therefore, 6(A/n) = iS UM 2 
A/n no — a(23°— t) 
and mo = 1.824 X 10~4, 
a = 493 <x 10", 
hence “an = — 0.27% per degree C. for air at 20°C. 
/n 
For oxygen 8(A/n) =— abst ’ 
rs A/n 1+ dt 
and a = 0.82, 
b = 3.665 <X 107°; 
hence HA/s) = — 0.33% per degree C. for oxygen at 20°C. 
/7 
For water vapor, 8(A/n) Pree. ae 
A/n no + a(t — 100) 
and a = 3.86 xX 107, 
“0 = 1.255 X 10™"; 
hence 8(A /n) = + 1.25% per degree C. for water vapor at 20°C. 





/n 
The fractional change in intercept is 
6T 


> ~ + 0.18% per degree at 20°C. 


(BVT) _ 
BJT 








411 


A SPECTROPHOTOMETRIC DETERMINATION OF EXHAUST 
GAS TEMPERATURES IN THE PULSE-JET ENGINE! 


By H. F. QuINN 


Abstract 


This paper describes a spectrophotometric method whereby instantaneous 
values of a variable flame temperature, in the particular case of nonluminous 
flames, may be determined and continuously recorded. 

This new technique, which depends upon the establishment of monochromatic 
black-body radiation conditions in the flame for a small region in the visible 
spectrum, involves the continuous measurement of radiation intensity in the 
above region, the intensity being, thereafter, correlated with the temperature 
of the flame. 

The problem of temperature measurement in the general case of nonluminous 
flames (flames which do not contain an appreciable amount of free carbon in the 
form of soot) is considered and a brief review of previous techniques employed 
for this purpose over the past 50 years is given. The basic theory and preliminary 
experimental justification of the present method are discussed. 

A description of the apparatus and the experimental arrangement used by 
the author in a specific application of the present method in the determination 
of the time variation of temperature in the exhaust flame of a pulse-jet engine 
is given. This includes details of a special type of spectrophotometer which 
employs a multiplier photocell as the radiation detecting and measuring element 
and, also, a “‘black-body”’ cavity constructed as a standard radiation source for 
the calibration of the former instrument. An original technique used to inves- 
tigate the emissivity of flames colored by alkali metal vapors is described and 
its application to the present problem shown. 

Finally, the measurable temperature range of the present apparatus is con- 
sidered together with the inherent limitations of the new method. 


Introduction 


Developments in the field of jet propulsion over the past 10 years have 
emphasized the necessity for the development of suitable techniques for the 
measurement and recording of instantaneous values of thrust, fuel flow, gas 
pressure, and temperature in typical jet propulsion motors. This process of 
instrumentation development has proceeded with great rapidity both in Great 
Britain and in the U.S.A. It is desired, in this paper, to report on one particular 
phase of the general development, namely the technique of measuring and 
recording continuously and rapidly varying temperatures in the exhaust gas 
stream of jet propulsion motors. To simplify the discussion, the paper has 
been divided into two parts, Part I of which deals with the basic physical 
principles and underlying theory of a new technique in flame temperature 
measurement, while Part II describes the application of the developed method 
in the case of a specific jet engine, the pulse-jet. 


The present problem, as originally assigned to the author (then engaged as 
a Research Associate of the University of New York) by the late Prof. J. K. L. 
MacDonald, Technical Director of a jet propulsion project (Project ‘‘Squid’’) 


1 Manuscript received in original form September 26, 1949, and, as revised, February 28, 
1950. 


Contribution from the Canadian Armament Research and Development Establishment, 
Valcartier, Que. 
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in the above University, may be summarized as follows:—‘‘It was desired to 
design and construct an instrument capable of providing an accurate record 
of varying flame temperatures at any specified point in a jet motor. This 
instrument was required to record these temperature variations up to a funda- 
mental frequency of 250 engine cycles per second and was to be fully responsive 
to harmonic frequency components up to the limit of 2500 cycles per second. 
The instrument was not to be subject to destructive action by the hot gases 
of the jet and was to be free from radiation or other sources of error. Finally, 
the insertion of any part of the instrument (if insertion were necessary) into 
the gas stream was not to lead to any additional perturbation of its mode of 
flow”. The above requirements were originally postulated with specific reference 
to one particular type of motor, the pulse-jet and, following some preliminary 
investigations, it was decided to limit the application of the instrument to the 
practically nonluminous gas stream at the exhaust end of one of the above 
engines. 


Part I 


A thorough search of the existing literature in the field of gas temperature 
measurements was made in the hope of finding an existing method of dealing 
with the problem or, at least, one which could be conveniently modified so to 
do. Of the mass of literature on the subject, dating from 1890 onwards, only 
one article gave promise of a method which might be employed after consider- 
able modification to meet the above stated requirements. A brief outline of the 
survey has been included here in so far as it is indicative of the many problems 
which are encountered in the temperature measurement of nonluminous 
flames. The very earliest methods of temperature measurement applicable to 
nonluminous flames required the introduction of thermocouples into the body 
of a stationary flame such as that of a Bunsen burner: many elaborate tech- 
niques were evolved to compensate for losses due both to radiation from the 
couple itself and to conduction along the leads. Typical experimental arrange- 
ments are described by Waggener (20), by Burkenbusch (2) and by Nicols (15). 
It is obvious that the time-lag due to the thermal capacity of any practical 
thermocouple element completely eliminated the use of such a device as far 
as the present problem was concerned. 


An original optical method of flame temperature determination was devised 
by Kurlbaum (9) for the luminous, carbonaceous type of flame but served as 
the basis for the method of spectral line reversal due to Fery (3) and inde- 
pendently to Kurlbaum and Schulze (10). An excellent description and theo- 
retical explanation of the latter method has been given by Lewis (11). While 
this method, which applies strictly to a nonluminous flame, has been used with 
considerable success by Awbery and Griffiths (5) and by Loomis and Periott 
(12) for the determination of static flame temperatures, it is essentially a 
static method in so far as it requires a visual matching of brightness by the 
experimenter and is, therefore, completely unsuitable for the determination of 
rapidly varying temperatures. It is noteworthy, however, that this method 
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has been extended in its application to the special case of nonluminous flames 
in an Otto cycle engine cylinder by the use of a stroboscopic arrangement 
whereby the brightness matching is effected at a specific point in the cycle 
of the engine; this is achieved by the use of a stroboscopic disk driven synch- 
ronously from the engine crankshaft, a set of holes cut into the disk allowing 
light from a comparison radiator to pass through the flame during a very short 
time interval at some specific point in the engine cycle. Procedures devised 
for this point-to-point temperature measurement are described by Hershey 
and Patton (6), and also by Rassweiler and Withrow (16). This modified line 
reversal technique is unfortunately inapplicable to the pulse-jet owing to the 
fact that not only are the individual cycles of the engine remarkably dissimilar 
but also because there is no rotating mechanical part of the engine with which to 
effect the synchronization of a stroboscopic disk. 


A further optical means of determining temperatures in the case of non- 
luminous flames depends upon the measurement of emission and absorption 
of radiation by the flame in the infrared region of the spectrum between the 
limits of 2 and 6 microns: this method was extensively investigated by Schmidt 
(18) for the particular case of the Bunsen flame. Schmidt was able to show 
conclusively that the Planck Radiation Law applied to certain molecular band 
radiation in the above region. In order to apply this method to the flame of the 
pulse-jet exhaust, it would first be necessary to examine the complete infrared 
spectrum from a gasoline flame burning in air under conditions which closely 
simulated those in the actual pulse-jet engine in the hope of finding some 
particular band radiation for which the above law applied: facilities for such a 
study were not available. No attempt was made, therefore, to explore this 
method any further. 


The experimental procedure finally adopted was first suggested by Mac- 
Donald following consideration of an original report by Graff (4) in which 
a description was given of another radiation method of temperature measure- 
ment applicable to the flames in an Otto cycle engine. Graff considered that, 
as the radiation properties of the intermediate and final products of combustion 
in the combustion chamber of such an engine were not sufficiently well defined 
and were, in part, unknown, it would be feasible to make exact temperature 
measurements using the radiation due to certain additive substances purposely 
introduced into the combustion chamber. Suitable additives were found to be 
the alkali metals introduced with the fuel in the form of organic salts reasonably 
soluble in gasoline. Graff quotes some preliminary investigations of Yosida (a 
reference which cannot be traced), wherein it was apparently shown that the 
above additives (when present in sufficient concentration in the flame) emit, 
within the limits of their resonance lines, radiation which is a function of 
temperature only, as given by the Planck formula: that is to say, monochro- 
matic black-body conditions could be established in the flame for the spectral 
regions of the resonance lines. It is at once apparent that, in general, the above 
condition is not attained, the radiation intensity due to the additive depending 
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not only upon the temperature but also upon the concentration of the additive 
in the particular portion of the flame examined. No mention is made by Graff 
of any investigations pertaining to the minimum amount of additive required 
to effect black-body conditions in the flames which he examined. 


The apparatus used by the above author in the application of his method to 
the Otto cycle engine flame consisted of a photoelectric cell used as a radiation 
detector, a color filter placed in front of the detector to eliminate band or con- 
tinuous spectral radiation from the flame, and an amplifier to amplify the 
output voltage of the photocell up to the level necessary to operate a cathode 
ray tube on the screen of which the intensity versus time record was presented. 
No mention is made of the calibration procedure employed or of the means of 
recording the deflection of the oscilloscope beam. In order to apply this 
method to the case of the pulse-jet exhaust flame, it was first necessary to ex- 
amine whether or not the concentration of the additive which it would be 
possible to inject into this engine could be increased beyond the minimum 
level required to provide monochromatic black-body conditions in the ex- 
haust flame. 


THEORETICAL 


The concept of temperature as applied to the particular case of the com- 
pletely or partially burned gases which comprise a flame is defined as corres- 
ponding to a particular state of complete statistical equilibrium in these gases. 
Such a state is described thus: collisions between the atoms or molecules and 
interaction with the radiation of the system, although changing continually 
the energy state of every individual particle, maintain, on the average, the 
randomness of direction of motion and the constancy of (a) the percentage of 
atoms or molecules possessing a specific velocity, (b) the quantity of molecules 
or atoms possessing a given quantum state of rotation, vibration and of elec- 
tronic excitation, and, finally, (c) the concentration of the dissociation pro- 
ducts. It is here assumed, following Lewis (11), that, in the present case of the 
exhaust flame, the equilibrium (a) is completely established owing to the 
rapidity of the molecular collision sequences. An additional factor in the 
statistical equilibrium must be postulated in the case of the present problem 
where an additive, sodium, is present in the exhaust gases, namely the con- 
stancy of the percentage of atoms of the additive possessing a given degree of 
electronic excitation. In this connection, it has been shown by Bauer (1) and 
by Awbery and Griffiths (5) that identical temperatures are obtained in prac- 
tice in the case of a specific flame by means of the line reversal technique irre- 
spective of the particular alkali metal used as the additive, that is, irrespective 
of the excitation potentials of the atoms. This experimental fact admits of no 
other interpretation than that a complete thermal equilibrium existed between 
the translational degrees of freedom of the gas molecules and the electronic 
degrees of freedom (excitation states) of the alkali metal atoms. Since the 
radiation intensities due to the atoms of the additive is intimately connected 
with the number of these atoms in the several excited states, it is possible to 
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associate a measured radiation intensity from a given concentration of these 
atoms with a specific temperature. In the case of the jet exhaust flame speci- 
fically considered here, this temperature will be defined to correspond to a 
state of statistical equilibrium between the translational degrees of freedom of 
the gas molecules and the electronic degree of freedom of the additive atoms. 
No attempt has been made to investigate the establishment of the (6) and (c) 
equilibria: such an attempt would be unnecessary in so far as modern gas 
dynamical calculations as applied to jet engines include the simplifying assump- 
tion that the gases are ‘“‘perfect’’. Hence, only a temperature corresponding to 
a statistical equilibrium in the translational degrees of freedom is of useful 
significance in such calculations. 


The intensity—temperature relation for the additive has been seen to alter, 
in general, with the concentration of the latter in the flame. Two limiting 
cases of the thermal excitation of atomic spectra may, however, be distin- 
guished to which the terms ‘Planck radiation” and “‘Boltzmann radiation” have 
been applied by Smit (19). Thus, consider a perfect black-body cavity, the walls 
of which are maintained at a temperature 7: let some means be provided 
whereby a controiled amount of some alkali additive (sodium, for example) can 
be admitted in the form of a gas. If now two small holes are bored axially into 
the cavity as shown in Fig. 1, it will be observed that the radiation emerging 


PLANCK RADIATION 


Resultant Line Emission 
(Weakened By Self-Absorption) 


Fic. 1. Theoretical black-body cavity containing a discrete line emitter. 


from the wall (direction 1 in the figure) will obey the Planck formula, whereas 
the radiation observed in the axial direction joining the two holes will give the 
line spectrum of the sodium gas, the intensities of the spectral lines depending 
upon the populations of the excitation levels according to the Boltzmann for- 
mula. This situation is represented in Fig. 2 where the nearly horizontal line 
represents the spectral energy distribution of the radiation proceeding from 
the wall, while the lower curve represents the energy distribution in the line 
spectrum observed along the axis. As the concentration of the additive is 
increased, so will the intensities of the several lines increase, being finally 
limited by self-absorption to the ‘‘ceiling’’ represented by the Planck curve: 
that is to say, monochromatic black-body conditions will have been established 
for the axial direction of observation for each of the spectral lines. 


It is possible, by means of a simple experiment, to show that monochromatic 
black-body conditions can likewise be established in the case of the “‘D”’ line 
resonance radiation from sodium atoms thermally excited in a flame. The 
experimental arrangement is as follows. A saturated solution of common salt is 
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Fic. 2. Distribution in intensity for Planck radiation from the cavity walls and Boltzmann 
radiation from the thermally excited line emitter. 
Planck Radiation 


2hv3 1 
I ia ar Fae ae 
”) = * Ge/KF oD) 
erg 


cm sec. sterad. sec.— 


aspirated under some 25 lb. air pressure into a Bunsen type burner of rectangular 
cross section together with ordinary coal gas used as fuel. A photocell used as 
a radiation detector receives ‘‘D”’ line radiation emanating from a 1 cm. square 
area of the flame at the front of the burner. Once steady burning of the flame 
has been established, an opaque barrier consisting of a sheet of “‘transite’’ is 
moved into several positions along the length of the burner, effectively reducing 
the thickness of the flame as seen by the detector as it moves towards the front 
of the burner. Fig. 3 shows a plot of typical results obtained in practice, 
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Fic. 3. Results of a preliminary test for monochromatic black-body conditions in a static flame. 
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intensity as measured by the detector being plotted against the distance of the 
barrier from the front of the burner. It is concluded from these results that, 
in the above case, monochromatic black-body conditions were established for 
the ‘‘D” line radiation in a 1 cm. thickness of the flame, the requisite amount 
of additive being obtained by the injection of the salt solution. Theoretically, 
the above conditions could be established in a lesser flame thickness but, in 
practice, it would be difficult to provide the necessary high concentration of 
additive. 


In any specific case of flame temperature measurement in terms of the 
radiation intensity of an alkali metal additive, it remains, therefore, to be 
determined whether or not a sufficient amount of the additive can be injected 
to establish a monochromatic black-body condition in the flame for the reson- 
ance radiation. This, in turn, presupposes some means whereby the establish- 
ment of the condition may be verified. If the condition can be shown to exist, 
then the intensity of the resonance radiation will follow the Planck law, and 
the calibration of a suitable spectrophotometer in terms of temperature against 
radiation intensity will readily be obtained by sighting the instrument into a 
standard black-body cavity. In particular, the case of the exhaust flame of a 
pulse-jet motor will now be considered, this flame, in the specific case of the 
small motor employed in these investigations, having a thickness of 4.14 cm. 
at the extreme end of the tail pipe. 


Part II 


APPLICATION OF THE SPECTROPHOTOMETRIC METHOD IN THE PARTICULAR 
CASE OF A PULSE-JET ENGINE 


Details of the operation of the pulsating jet engine, together with a com- 
prehensive review of the history of its development, have been described by 
Miller (14); consequently a description of the engine and its properties is 
omitted here. According to the above author, satisfactory instrumentation 
techniques for measuring the time average of the thrust, fuel flow, and body 
temperature of the engine were developed after minor modification of existing 
methods but, as of the summer of 1946, development of satisfactory instru- 
ments to measure and to record instantaneous values of thrust, pressure, and 
temperature had proved to be a slow and difficult process primarily owing to 
the high operating rate of the engine (up to some 250 complete engine cycles 
per second). 


Great impetus was given to the development of adequate instrumentation 
not only by the pressing requirements for such devices in the empirical de- 
velopment of these engines at several laboratories in the U.S.A. but also by 
the need for actual numerical values of instantaneous pressure, temperature, 
gas density, and gas velocity for substitution in the theoretically derived hydro- 
thermodynamical equations for the pulse-jet which had been developed by 
MacDonald and his school (13) in their one dimensional mathematical treat- 
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ment of pulse-jet operation. At this time, therefore, among the other instru- 
mentation requirements, it became imperative to develop some method of 
measuring and recording the rapid cyclic temperature variations at several 
points in the main combustion chamber, the transition cone, the tail pipe, and 
finally in the exhaust flame itself. 


Information received from the Naval Research Station at Annapolis, where 
continuous performance tests of pulse-jet engines were being carried out at 
this time, indicated that the flames in all portions of the particular types tested 
there were nonluminous (as opposed to the luminous, soot-containing type of 
flame). Accordingly, the instrument requirements were initially limited to a 
device which would meet the previously stated specification in the case of 
nonluminous flames only. It was not until the construction and testing of a 
glass-walled pulse jet had finally been accomplished at New York that the 
flames both in the combustion chamber and in the tail pipe were revealed to 
be brightly luminous, only the exhaust flame being nonluminous (a faint bluish 
flame). At this point, the research program diverged, a group at New York 
commenced investigation of methods of temperature measurement for luminous 
flames while the author, having at this time returned to McGill University, 
continued the original research as applied to the case of the nonluminous 


exhaust flames. 
METHOD 


The general arrangement of the apparatus as employed in the pulse-jet 
application is illustrated in Fig. 4. The sodium additive was, at all times, in- 
jected with the air stream through the intake valves of the motor in the form 
of a finely powdered inorganic salt by the use of a small blower located some 
distance away from the intake valve bank. This method of additive injection 
was considered far superior to the water-spray injection technique used by 
other investigators, as there is less probability of the operation of the engine 
being adversely affected by the addition of a dry powder than by the consid- 
erable amount of water required to effect adequate aqueous injection of the 
additive; cooling of the flame due to the absorption of the heat of vaporization 
of the water is also eliminated. 


The arrangement of the optical and electrical components of the simple type 
of spectrophotometer employed is shown in Figs. 5 and 6 respectively. An 
image of the source (a small hole of 13/64 in. diameter pierced into the blast 
and radiation screen directly opposite the extremity of the tail pipe of the 
motor) is formed by an auxiliary biconvex lens on the spectrometer slit. The 
prism P is set in the position of minimum deviation for the “‘D” line image 
which is formed by the lens Ze at the knife-edge jaws of a second adjustable 
slit, Ao. The latter opening is adjusted until the bright image of the ‘‘D”’ lines 
just disappears through it, the light in the spectral region so selected being 
received by the lens L; whose focal plane is at the second slit. Finally, the 
light emerges to uniformly flood the photo-cathode of the detector element, a 
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Fic. 4. General arrangement of the apparatus in the pulse-jet application. 


To Te 


Fic. 5. Optical system of the spectrophotometer. 


L,— Collimator lens; focal length, 7.525 in. L2— Telescope lens; focal length, 7.525 in. 
L3— Objective lens; focal length, 1.746 in. A,— Spectrometer slit. 
Az— Selector slit. P — 60° crown glass prism. 
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Fic. 6. Electrical circuit of the spectrophotometer. 


R,— 1.05 M ohms T:— To microammeter 

R2 to Rio— 100 K ohms T:— To d-c. amplifier 

Ru— 50 K ohms S:— Change-over switch. 

Rizx— 200 K ohms E — Supply voltage, 2000 v. d-c. regulated. 


standard electron multiplier photocell (R.C.A. type 931-A). The requisite 
voltages for the several dynodes of the tube are obtained as shown from a 
potentiometer, R; to Ru, placed across the output of a 2000 v. regulated d-c. 
power supply, £. The entire spectrophotometer is housed in a light-tight metal 
box, the latter also protecting the detector from the considerable blast and 
vibration from the jet engine. Further protection from these effects is also 
obtained by mounting the instrument on a separate heavy table. 


A check on both the optical filtering efficiency of the above system and also 
on the linearity of the photocell was obtained by the use of a calibrated broad- 
filament tungsten lamp as a light source: the phototube currents corresponding 
to a series of different lamp currents were recorded. The tungsten lamp had 
been previously calibrated in another experiment and the relation between its 
apparent black-body temperature and heating current was known. Provided 
the filtering was adequate and the tube response linear, a plot of the natural 
logarithm of fhe measured intensity in the immediate spectral region of the 
“D” lines against the reciprocal of the absolute temperature should yield a 
straight line, that is, the Wien law equation should hold. This plot is shown in 
Fig. 7: taking the ‘‘mean wave length”’ of 5893 Angstrom units and the accepted 
values for the Planck and Boltzmann constants, the theoretical slope of the 
straight line was computed, giving the figure —2.48. The experimental value 
obtained from the above graph was —2.57, leading to an error of just under 
4% which, considering the roughness of the measurements, was taken to con- 
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firm both a satisfactory filtering of the incident light and a satisfactory linear 
response of the tube. 


Light Source- Calibrated Tungsten Filament 


~ ol 


Log I (micro-amperes in 931-A Load ) 


5 
1/T degrees K"x 104 
Fic. 7. Results of a preliminary test of the spectrophotometer (application of the Wien law 
in the form— — 


log, I = log, Ci -& 7° 


A single stage of push-pull d—c. amplification was used between the photo- 
cell and the cathode-ray oscilloscope. It was experimentally determined that, 
by increasing the level of the light incident upon the multiplier tube and 
by the use of a 200,000 ohm resistor as a load, it was possible to obtain an 
output voltage adequate to operate the push-pull single stage d-c. amplifier 
which is employed as a fourth stage in a Dumont type 208-B oscilloscope. This 
amplifier stage has a frequency response flat from 0 to 100,000 c.p.s. and has 
its output terminals directly connected to the deflection plates of the cathode 
ray tube. Frequent tests demonstrated the excellent stability of this system, 
and its operation was completely adequate in all respects. 


A procedure for determining the existence of monochromatic black-body 
conditions in the case of a static, steadily burning flame has already been 
described: in so far as this method requires the insertion of a barrier into the 
flame, it is inapplicable to the jet exhaust flame—nothing must be inserted into 
the latter which might perturb the normal mode of flow of the exhaust gases. 
To overcome this difficulty, it was decided to assume the applicability of the 
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Kirchhoff law for the ‘‘D’’ line radiation (17) and to test the flame for unity 
monochromatic absorptivity. This was arranged as follows. Radiation from a 
small area of the flame diametrically opposite the source hole in the screen was 
received by a small plane mirror and reflected back into the flame towards the 
spectrophotometer. If the flame possessed unity absorptivity for the ‘‘D’’ line 
radiation, it was considered that identical intensity values would be recorded 
by the photometer whether the mirror was exposed to the radiation or not. 
Preliminary tests on static flames demonstrated the effectiveness of the method. 


In the application of this test to the exhaust flame, it was decided to provide 
a satisfactory basis for comparison of intensities measured both with the mirror 
exposed and screened in the following manner. A rotating metal disk having 
four open sectors cut into it was arranged to spin in front of the mirror, thereby 
alternately exposing and screening the latter during a period corresponding to 
four complete engine cycles. It was, thereby, anticipated that all cycles would 
show approximately the same maximum amplitude of intensity and general 
intensity—time relation if unity absorptivity were established, while alternate 
sets of four cycles would show marked differences if this condition were not 
attained. A provision was made whereby the exact period during which the 
mirror was screened could be indicated on the final intensity—-time record by 
the use of a cam fitted to the drive shaft of the disk: this cam was arranged to 


Fic. 8. Emissivity test arrangement. 
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close a sensitive microswitch whenever a segment of the disk completely 
screened the mirror. The switch itself was connected to the external contactor 
terminals of a General Radio “‘Strobotac” which provided a brilliant flash of 
light each time complete screening of the mirror was accomplished. These light 
flashes fell on the face of the cathode ray tube screen and were recorded simul- 
taneously with the deflection of the beam. The sector disk, the driving motor 
(operated at 3150 r.p.m.) and the plane mirror are fully illustrated in Fig. 8. 


The calibration method employed a cylindrical black-body cavity heated 
electrically to incandescence by a concentric graphite heating resistor as a 
standard radiation source. The complete unit is shown in Fig. 9. The cavity 


Fic. 9—A 


Fics. 9-A AND 9-B. The black-body cavity and furnace. 

Fig. 9-A. 

W — Water cooling pipes E,— High-current clamp electrodes 

Y — Carbon end sections C — Black-body cavity 

G — Graphite tube Z — Zirconium oxide tube 

B — Cavity end insert D — Radiation breaker 

V — Variable transformer 
Fig. 9-B. 

M — 0-15 amp. a-c. meter 

X — High-current transformer 

Re — Resistance (Tube G) 
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itself was constructed by core-drilling out of a block of solid zirconium oxide, 
this refractory material having a melting point in the vicinity of 2500°C..A 
central pivot arrangement together with the provision of flexible leads allowed 
the entire furnace to be swung from one fixed position where it could be viewed 
by the spectrophotometer through an angle of 30 degrees into another position 
where it could be viewed by a Leeds and Northrup type 8622 optical pyro- 
meter. Thus a measured radiation intensity could be rapidly correlated with 
a specific black-body temperature of the cavity. 


The circular aperture of the cavity was chosen to be identical with the source 
hole in the blast screen; hence the images formed on the spectrometer slit both 
from the flame itself and from the cavity were of equal size. It is fully realized 
that perfect black-body conditions were not attained in the cavity owing to 
the nonuniform electrical heating system employed and to the extremely poor 
thermal conductivity of the zirconium oxide which introduced large temper- 
ature gradients in the walls. In the absence of a tungsten cavity, however, the 
above arrangement was the most satisfactory; apart from the necessity for 
frequent replacement of the heating element the entire unit gave trouble-free 
service throughout its period of use. 


Recording of the intensity-time variations was achieved throughout the 
entire period of research by the use of the conventional cathode ray oscillo- 
scope — continuous film camera arrangement, the output signal of the amplifier 
being presented in the form of an ‘‘X”’ deflection of the electron beam. The 


oscilloscope was a standard Dumont type 208-B model, the camera a General 
Radio type 651-A-E incorporating an f/1.5 lens. Timing marks were impressed 
on the film by the Strobotac arrangement as previously described. 


The procedure adopted during a single test run on the engine was as follows. 
A period of one-half hour was allowed in order for the electronic circuits to 
become stabilized: during this waiting period, the camera was loaded with 
50 ft. of film, the a-c. mains voltage checked and recorded, the main optical 
chain of the spectrophotometer aligned and the additive injection blower loaded 
with finely powdered sodium chloride and set in position facing the intake 
valves of the jet. Owing to the extreme sensitivity of the multiplier photocell, 
all the tests were run in complete darkness, apart from the small, flashing light 
of the Strobotac. As soon as the room was darkened, the high voltage supply 
to the cell was switched on, together with the ‘‘Beam”’ switch of the oscillo- 
scope, and the disk motor was started. As the entire film length was run through 
the camera in approximately 1 1/2 sec., it was imperative that the switching 
on of the camera be effected by one operator who carefully observed the de- 
flection on the oscilloscope screen while another operator assumed respon- 
sibility for both the starting of the engine and the injection of the additive. 
At the start of the test run, the latter operator switched on the starting spark 
and aspirated fuel into the engine combustion chamber, using an auxiliary 
starting compressed air supply. Once stable operation of the jet had been 
attained over a five-second period of time, the jet operator threw on the blower- 
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injector switch, forcing a stream of salt dust into the valve bank of the engine. 
Meanwhile the camera operator observed the oscilloscope screen and, when 
the single spot on the latter expanded into a trace some three inches long 
consequent upon the injection of the additive, made the main switch to the 
camera motors. The entire procedure from the initial starting of the engine to 
the final switching off of the camera occupied just under 15 sec. 


Completion of this recording was followed by a calibration of the photometer 
as follows. The engine, the emissivity test apparatus, and the blast screen were 
removed and the furnace placed in position. The latter was slowly raised 
through a series of temperatures from 800 to 1500°C. (these temperatures being 
measured by the optical pyrometer) and the deflections corresponding to 
specific cavity temperatures were recorded. The results of a typical calibration 
run appear in Table I. This data enabled the curve of Fig. 10 to be con- 
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Light Source - Black Body Cavity 


Slope (Calculated) - 2.46 
Slope (Measured ) - 2.43 


O- Denotes Run No. ! 


e-Denotes Run No. 2 


Log, C.R.T. Beam Deflection 


2 ot @ 
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Spectrophotometer calibration curves (application of the Wien law). 
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TABLE I 





Cavity temperature 
oe Oscilloscope deflection, in. 
°F °G °K 


2300 1262 1535 
2262 1240 1513 
2242 1230 1503 
2230 1223 1496 
2220 1219 1492 
2110 1157 1430 
2020 1105 1378 
1965 1072 1345 
1880 1028 1301 
1790 976 1249 
1760 959 1232 
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Fic. 11. Corrected calibration curve for the recording oscilloscope. 
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structed where it is to be particularly noted that the measured slope of the 
curve differs by only 1.2% from the theoretically calculated value, thus indi- 
cating the satisfactory achievement of black-body conditions in the cavity. 
The final calibration curve for the spectrophotometer-—oscilloscope combination 
appears in Fig. 11. It will be noted that the thermally excited ‘‘D” line radia- 
tion, as detected by the present apparatus, tends to a practical cutoff in the 
vicinity of 950°C. 


The conversion from the intensity-time record on the film to the final 
temperature-time results was effected as follows. Two marks were provided 
on the face of the cathode ray tube, the distance between these being carefully 
measured, using dividers and a steel rule. Finally the two marks were photo- 
graphed and the negatives obtained placed in a Kodak Precision Enlarger 
which was then adjusted until the distance between the images of the two 
marks projected on a plane white surface equalled the distance originally mea- 
sured on the face of the tube. Having made this adjustment of the enlarger, 
it was locked in this position: thereafter it was simply necessary to project 
any particular film record onto the white surface and, by the use of the cali- 
bration curve of Fig. 11 to read off the temperature-time results. 


Results 


The following results have been derived from a study of some 600 ft. of film 
constituting the intensity-time records for 12 separate test runs of the jet 
engine. First, it was immediately apparent on examining the records that. a 
general dissimilarity existed in the intensity-time wave forms from one cycle 
of the engine to another. This is exemplified in Fig. 12 where several cycles are 
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Fic. 12. Representative set of successive intensity-time cycles obtained in the pulse-jet appli- 
cation. 


shown, these being recorded before the camera reached its maximum operating 
speed of 25 f.p.s. Indeed, when all the factors which influence the operation of 
this type of engine are reviewed particularly the nature of the fuel injection 
system, it was considered remarkable that any consistency among the wave 
forms would be obtained: some definite consistency has, however, been defin- 
itely established in this research. 


In examining the film records, a search was made for two particular points 
of similarity in the cycles: (1) equality of the intensity maxima and, (2) simil- 
arity of the general configuration of the wave forms. It was discovered that 
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one special configuration occurred at irregular intervals throughout the records: 
this configuration is illustrated in the rough sketch below. 


Intensity (units arbitrary ) 





Time (units arbitrary) 


It is to be noted that the test for the establishment of black-body conditions 
in the exhaust flame did not give positive results in that no alternate sets of 
four large and four small intensity-time wave forms were ever observed, the 
maximum intensity and general configuration varying widely from cycle to 
cycle. It is interesting to note, however, that the records show that the above- 
mentioned scattered similar cycles were recorded at times when the mirror was 
screened and also when it was exposed. 


The existence of these recurrent similar cyles admits of two separate inter- 
pretations: (1) there is a fairly definite and consistent temperature-time varia- 
tion in the exhaust gases and, in the case of the above cycles, black-body 
conditions were established in the flame or, (2) the temperature-time relation 
varies widely from cycle to cycle, black-body conditions were never established 
at any time in the flame and the above similarity results only from fortuitous 
combination of several factors. These factors are (a) the particular temperature— 
time variation for a given cycle, (b) the particular concentrations of the addi- 
tive present in the exhaust flame at all times during the cycle and, (c) the effect 
of additional light received by the spectrophotometer in those cases where the 
mirror was exposed. The former interpretation, being the most probable, has 
been accepted here. 


While an interpretation of the consistent cycles has been given, the wide 
dissimilarity observed from cycle to cycle remains to be explained. Studies of 
the behavior of the engine, particularly in the cases where an especially finely 
powdered additive was injected, indicated that the normal combustion processes 
in the engine could be adversely affected by the addition of an excess of salt 
powder: such a disturbance would explain the widespread variations from one 
cycle to the next. It is considered, however, that at some particular intervals 
during the recording periods, the concentration of the additive present in the 
engine was adequate to effect the requisite conditions in the exhaust flame 
without seriously perturbing the main combustion process and that the above 
consistent cycles were recorded at such times. 
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The return of the intensity-time wave forms to the zero level over a large 
part of the operating cycle is due to two separate factors. The present apparatus 
is unable to detect thermally excited ‘‘D”’ line radiation emitted by a source 
whose temperature is less than 950°C. In the absence of evidence from other 
investigations revealing a more surprising reason, it might be reasonably as- 
sumed that, at the periods of zero intensity, the exhaust gases simply fell to 
temperatures below the above limit. It has, however, been definitely established 
by Hudson and Hett (7) in their investigations of the pulse-jet exhaust flame 
by means of high-speed photographic techniques that, following the initial 
outflow of the exhaust gases shortly after the beginning of a particular engine 
cycle, a period of flow reversal occurs wherein the exhaust flame appears to 
contract in diameter and then to vanish backwards into the tail pipe of the 
engine—there is a reversal of mass-flow. It is considered that both of the above 
effects operated to give the wave forms observed in this research; the exhaust 
gases probably fell first to a much lower temperature than 950°C. and then 
were effectively swept out of the field of view of the spectrophotometer, being 
replaced by atmospheric air at room temperature. 


MIRROR EXPOSED MIRROR SCREENED 


MIRROR SCREENED MIRROR EXPOSED 


MIRROR SCREENED F MIRROR EXPOSED 


Selected similar cycles obtained from records in the pulse-jet application. 
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While the lower limit of temperature which the present apparatus can 
measure has been seen to be 950°C., there is no theoretical upper limit to the 
measurable temperature range in so far as the maximum temperature which 
can be recorded is controlled only by the saturation current of the particular 
photocell employed, the ability of the amplifier to handle a large input voltage 
signal and, finally, the size of the cathode ray tube required to present the 
deflection corresponding to the maximum temperature. The apparatus here 
was limited in practice to an upper temperature limit of 1350°C. 


The appearance of the ‘‘secondary maxima”’ near the start of the intensity— 
time curves was initially regarded as being extremely peculiar in so far as it 
was difficult to understand why a sudden drop in temperature of the exhaust 
gases should occur shortly after the initial outflow. This peculiar phenomenon 
has, however, been strikingly correlated with some early high-speed photo- 
graphic results of Hudson and Hett (8). These investigators showed that the 
flame propagation in the pulse-jet takes place in the following manner:—upon 
ignition of the air-fuel mixture in the combustion chamber, a small knot of 
flame (termed a ‘‘flambeau’”’ by Hett) is rapidly propagated throughout the 
length of the engine, finally emerging at the end of the tail pipe; this preliminary 
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Fic. 14. Temperature-time relations for the particular cycles 13-A and 13-E. 
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“‘flambeau”’ is followed after an interval of no visible burning, by a marked 
bank of flame traveling some four inches behind the flambeau. The occurrence 
of the above secondary maximum is, therefore, interpreted to correspond to 
the arrival at the extremity of the tail pipe of the high temperature section of 
the flambeau, while the succeeding upward rise to the main maximum is taken 
to correspond to the arrival of the main body of flame. 


When the several consistent cycles which are shown in Figs. 13-A to 13-F 
are considered, the logarithmic relation between the recorded intensity and 
the actual temperature must be borne in mind: in particular, it must be realized 
that, at the higher temperatures, a small change in temperature results in a 
very disproportionate change in intensity. To further illustrate this point, the 
temperature-time graphs obtained from the particular cycles of Figs. 13-A 
and 13-E are shown in Fig. 14. Both curves indicate a secondary maximum 
occurring 1.0 msec. after the beginning of the initial upward temperature 
gradient; moreover, the general configuration of the curves and their maximum 
temperatures are strikingly similar. Fig. 15 shows similar temperature-time 
curves for the cycles of Figs. 13-B-and 13-C. As minor differences still exist 
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between the above curves, it has been decided to select Curve 1 of Fig. 14 as 
representing an average relation, and this particular curve has been taken to 
represent the considered temperature-time relation for the exhaust flame. 
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FAST NEUTRONS FROM THE T + D AND T + Li REACTIONS! 







By E. ALMQvVIST 


Abstract 


The fast neutrons produced by the triton bombardment of LiF, LisCOs, 
LiOH.H:0, LiOD, and LiOD.D20 + D.O have been detected, and the thick 
target yields and their slowing down lengths in graphite measured. The tritons 
were produced by irradiating the compounds in the high thermal neutron flux 
of the thermal column of the NRX pile. The spatial distributions of the fast 
neutrons from the T + D and T + Li reactions were measured in the graphite 
with indium resonance detectors. The numbers of fast neutrons emitted per 
second were obtained by comparing the space integrals of the indium activity 
with that measured when a standard Ra-a-Be source of neutrons was in the 
position of the lithium compounds. The number of incident tritons was estimated 
from the measured flux of thermal neutrons at the lithium compound and the 
total capture cross section of the lithium. The thick target yields were thereby 
derived. The average cross sections for neutron production by 2.6 Mev. tritons 
are 1.5 X 10-*‘cm.? per atom of lithium and 0.41 X 10-*‘cm.? per atom of deu- 
terium. The slowing down lengths in graphite of density 1.64 gm. per cc. are 
19.2 and 23.4 cm. for the neutrons from the lithium and deuterium respectively. 
The smallness of the latter indicates that inelastic scattering of these 14 Mev. 
neutrons occurs in graphite. 

























Introduction 










The bombardment of light nuclei with tritons is expected to yield very 
energetic neutrons owing to the large positive energy balance of (T, ) reactions 
with most light elements. The reaction T(d, a)” has been studied at low energy 
(5) and the excitation function and cross section measured up to 2.5 Mev. 
deuteron energy (4, 17) by detecting the a-particles. Even at low bombarding 
energies this reaction has a prolific yield, and it is expected that triton reactions 
with other light elements will be useful sources of fast neutrons when tritium 
becomes available in sufficient quantities to be used in accelerators. 









Another means of obtaining fast tritons is offered by the high neutron flux 
in the NRX pile. The tritons are released in the reaction Li®(m, a)T for which 
the cross section for thermal neutrons is 870 barns. It has been possible to 
study the secondary reactions due to the tritons producing fast neutrons in 
the graphite thermal column of the pile. Here the naturally occurring level of 
fast neutrons is very low. 









This report describes an experiment to measure the thick target yields and 
slowing down lengths in graphite of the fast neutrons produced by triton 
bombardment of lithium compounds, some containing deuterium. In addition 
to the fast neutrons from the T + D reaction, neutrons from a T + Li reaction 
were detected and studied. The latter appears not to have been reported pre- 
viously in the literature. 










When a lithium compound is irradiated with thermal neutrons, tritons of 





1 Manuscript received March 20, 1950. : 
Contribution from the Nuclear Physics Branch, Atomic Energy Project, National Research 
Council of Canada, Chalk River, Ontario, Canada. Issued as N.R.C. No. 2162. 
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2.65 Mev. energy and a-particles of 2.00 Mev. energy are produced by the 
reaction Li(n, a)T (6, 7). These particles bombard the nuclei of the lithium 
compound and may produce reactions of the types: 


Th ee eee. (a) 


and 
a + Z4& —————> (Z + 2)4 +? + 03+ Q, (b) 


which give fast neutrons if the energy balance Q is positive. For the common 
isotopes of the light elements, except boron and beryllium, reactions of type (d) 
are endothermic (see Table II) and will not occur with a-particles of 2.0 Mev. 
energy, or if they are energetically possible the yield is small (13) owing to the 
short range of the a-particle in the material of the target and the effect of the 
potential barrier of the nucleus. Therefore, the fast neutrons produced can be 
assigned to triton reactions of type (a), and the data obtained give the thick 
target yields for 2.6 Mev. tritons on the compounds used. By using a detector 
which is insensitive to slow neutrons, these fast neutrons can be detected and 
information about their yield and energy obtained in the presence of a high 
thermal neutron flux. 


In this experiment calibrated indium resonance detectors were used to detect 
the fast neutrons after they had been slowed down to indium resonance energy, 
and the number of fast neutrons emitted per second by the source obtained 
by graphical evaluation of the space integral of the indium activity (1, 2). 
The number of incident tritons was calculated from the thermal neutron flux 
at the surface of the source, the amount of lithium in the source, and the 
thermal neutron capture cross section of lithium. 


Experimental Arrangement 


The thermal column of the Chalk River pile was used as a large block of 
scattering material with a high flux of thermal neutrons, as shown in Fig. 1. 
The lithium compound was loaded into a spherical shell container, 12 cm. in 
outside diameter, which was placed on the central line of the column 152 cm. 
from the front face. Each container consisted of two concentric shells made of 
1/32 in. stainless steel with 1 cm. space between them to hold the lithium 
compound, and a solid graphite sphere to fill the central space. With the pile 
operating, the thermal neutron flux at the surface of the shell was measured 
with four calibrated manganese foils fastened with cellulose tape to the 
graphite blocks in contact with the container. The thermal neutron flux was 
also monitored at the front face of the thermal column by indium and man- 
ganese detectors and the pile power monitored at an experimental hole so 
that separate runs could be normalized. The central string of graphite blocks 
were slotted to hold the indium resonance detectors, which were indium foils 
3.3 by 5.7 cm. in area and 113 mgm. per cm.? thick contained in disk-shaped 
cadmium boxes. These were made in two halves each that screwed together so 
as to have a wall thickness 2 mm. of cadmium on all sides. The detectors were 
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calibrated by putting a Ra-a-Be source of known neutron emission in place of 
the shell during a shutdown period of the pile and obtaining the space integral 
of the activity as described below. 
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Fic. 1. Experimental arrangement in thermal column of NRX pile. 











To minimize the effect of one indium detector on the adjacent ones they 
were loaded in alternate slots for each run, taking care always to place them 
with the same side facing the source and filling the remaining slots with graphite 
disks. After the activation these detectors were wrapped around a thin-walled 
aluminum £-counter and counted, first with one side facing the counter and 
then the other to minimize the orientation effect caused by the gradient of the 
neutron density. The background counting rate was subtracted and the read- 
ings corrected to zero decay and to saturation irradiation, taking 54.5 min. as 
the half-period of the indium activity. To calibrate the counters for sensitivity, 
six identical indium foils were irradiated in fixed positions in a paraffin wax 
box containing a 1.2 gm. Ra-a-Be source of neutrons and their activities 
measured on the counters each day. The counting losses at high counting rates 
were estimated from the measured decay of actinium active deposit, the half- 
period of which is known accurately as 36.1 min. (14), and corrections applied 
where necessary. 





















Calibration of Indium Resonance Detectors 


If a point source of fast neutrons is surrounded by a large volume of a slowing 
down medium in which there is no capture of neutrons above thermal energy, 









then co 
; Q, = sr q(E, r)r'dr, 
0 
where Q, = emission from the point source in neutrons per second, 
r = radial distance from the source, 
and g(£, r) = number of neutrons per second passing through the energy 






E per cc. at distance r from the source. 
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If a detector, sensitive only to neutrons of essentially one energy Eo is used, 
e.g., an indium resonance detector, then its saturation activity A(r) is pro- 
portional to the number of neutrons per second that pass through the energy Eo 
per cc. at the detector: 


i.€., q(Eo, r) = RA(r). 


Therefore Q, = 4xk [ A(r)r’dr = ko | A(r)rdr, 
J0 0 
where ky = 4k, 


and A(r) = activity of the detector placed at a distance 7 from the source. 


For any given source, A(r) can be measured by experiment at various values 
of r and the integral evaluated graphically in arbitrary units. To obtain the 
constant kp it is necessary to use a source of known neutron emission Q, and 
of sufficient strength to give the required statistical accuracy in the counting 
rates A(r) at large values of r. Once ko is known, the foils can be used to measure 
the emission of fast neutrons from another source. 


Two sources for which the absolute neutron emission had been determined 
were available. Source No. 1 consisted of 600 mgm. of radium surrounded by a 
beryllium capsule weighing about 8 gm. By means of the spatial integration 
method (1, 2) in water the neutron emission of this source had been compared 
with that of the same radium surrounded by a similar beryllium capsule. 
Gluckaiif (9) by a method developed by Paneth and Gluckaiif (10) then 
measured the amount of helium produced by the (y, 2) reaction on the beryl- 
lium in the latter capsule. This measurement, combined with the comparison 
of the emissions from the two beryllium capsules, gave an emission of 3.03 X 
10° neutrons per sec. from source No. 1 with an error estimated at + 3%. 
Source No. 2 was a 200 mgm. Ra-a-Be source that had been calibrated in terms 
of the standard source at the Argonne National Laboratory, Chicago. The 
absolute neutron emission of the standard source had been obtained by Seidl 
and Harris (16). The emission of source No. 2 was given as 2.8 X 10° neutrons 
per second, + 7%. 


As these standards were too weak to give good counting rates at large dis- 
tances from the source in the graphite, a 1.2 gm. Ra-a-Be source was calibrated 
in terms of the standards by comparing the activities induced in a tank of 
manganese sulphate solution (3, 19) irradiated by each of the three sources in 
turn. The tank was 60 cm. in diameter and filled to a level of 60 cm. depth with 
a nearly saturated solution of manganese sulphate (density of solution = 1.37 
gm. per cc.). Each neutron source was suspended by a string in the center of 
the liquid overnight. After removing the source from the tank the manganese 
y-ray activity was measured by observing the rate of deflection of a quartz- 
fiber electroscope connected to a large ion chamber filled to a pressure of 20 atm. 
with argon and immersed in the solution. A correction of 5% for loss of neutrons 
through the wall of the tank was applied in the case of the Ra-a-Be sources. 
Owing to the shorter slowing down length of the neutrons from the Ra-y-Be 
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source [2.5 cm. in water as compared with 6.8 cm. for Ra-a-Be neutrons (12)] 
no correction to the readings obtained with the photo-neutron source was 
required. 







The results obtained were 

1.2 gm. source = 64(+ 4%) X No. 1 = 64 X 3.03 XK 10°= 19.4 X 108 neu- 
trons per sec. 

1.2 gm. source = 6.42(+ 2%) X No. 2 = 6.42 X 2.8 X 10°= 18.0 X 10° neu- 
trons per sec. 

Mean value of 1.2 gm. Ra-a-Be emission = 18.7 X 10° neutrons per sec. 















Using the 1.2 gm. neutron source in the chosen position in the thermal 
column with the pile shut down, the value 





ce 
| A(r)r'dr = 2.35 X 10° counts per min.: X cm.’ was obtained. 
0 





Hence the detector constant 


18.7 X 10° ~.,..,-». neutrons per second 
9 = ———— = 7.96 X 10° * —_— ctl orate = 
2.35 X 108 counts per minute X cm.’ 















Calibration of Manganese Monitors of Thermal Neutron Flux 






There were available manganese foils, 1 cm. square, that had been calibrated 
by irradiating them in a thermal neutron flux the absolute value of which had 
been determined by means of a small, BF;-filled, pulse ion chamber of known 
counting volume and boron content (8). The activity of the foils was measured 
by counting with a Geiger—Miiller counter in standard geometry, and cali- 
brated for sensitivity by a standard uranium foil. 










The calibrated foils were not suitable for measuring the strong neutron 
fluxes encountered in the thermal column, as they became too active for the 
available measuring equipment. Instead, small disks of manganese foil, 1/8 in. 
in diameter, were stamped out and pressed into small | cm. squares of 3/64 in. 
aluminum sheet, which could be mounted in the standard foil holders. The 
small foils were calibrated for flux measurement by irradiating them simul- 
taneously with the standard foils in a paraffin box with a 1.2 gm. Ra-a-Be 
source of neutrons and counting them on the standard counters. This was 
repeated until every foil had been irradiated in each position in the box, and 
the mean counting rate was associated with the average flux of thermal neu- 
trons. This comparison showed a variation of + 4% in sensitivity among the 


small foils. 


When used in the thermal column the small foils could be measured con- 
veniently on a TQQ electroscope* but were too active to count in the standard 
geometry, and it was therefore necessary to obtain the conversion factor to 
change rate of deflection in divisions per minute on the electroscope to counts 
per minute on a standard counter. This was done by following the decay of a 



















* This type of electroscope was designed by Dr. H. Carmichael of this laboratory. 
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foil initially on the electroscope, then on an ion chamber connected to a 
vibrating reed electrometer, and when the activity had decayed sufficiently 
to count it on three of the standard counters. It was found that 1 div. per min. 
on the electroscope = 186 mv. per min. on the ion chamber and that 1 mv. 
per min. on the ion chamber = 259.5 counts per min. A blank aluminum 
mount was irradiated and counted simultaneously and its activity subtracted 
to get the activity of the manganese without a holder. The readings, converted 
to divisions per minute on the electroscope, are plotted on a logarithmic scale 
against time in Fig. 2. The slope yields a half-period of 2.60 + 0.01 hr., in good 
agreement with the accepted value of 2.59 hr. (15) for Mn®. 
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Fic. 2. Decay of activity of manganese foils used to monitor thermal neutron flux. Half-period 
= 2.60 + 0.01 hr. 


From the above intercomparisons the small manganese foils were calibrated 
to measure high thermal neutron fluxes with the aid of a TQQ electroscope. 


Results 


The following lithium compounds were irradiated: lithium fluoride, lithium 
carbonate, lithium hydroxide, lithium deuteroxide paste (LiOD.D.O + D.O), 
and lithium deuteroxide.* In order to determine the background due to epi- 
cadmium neutrons from the pile, runs were made with a mixture of boric acid 
and graphite, adjusted to give about the same total capture cross section as 
the lithium compounds (all approximately 300 cm.*). To check the yield of 
photo-neutrons from the deuterium in the deuteroxide a run was made with 
Ca(OD):2. The counting rates normalized as described in the following para- 
graph are plotted against the distance from the center of the source in Fig. 3. 

* The lithium fluoride and lithium carbonate were N. F. grade chemicals supplied by Mallin- 
ckrodt Co. and the lithium hydroxide C.P. grade. The lithium deuteroxide was made by the Chemical 
Control Branch and analysed as 84.4% LiOD, 5% moisture (mostly D2O), 2.8% CaCO; and 8% 


insoluble. The LiOD.D:O + D2O was obtained by adding D,0 to this. Analysis showed the LiOD 
to be 72 atom %D. The D2O added to form the LiOD.D.O + D:0 was better than 98 atom %D. 
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Fic. 3. Detector activity, A(r), vs. distance, r, from the source. The curve for boric acid 
represents the activity due to the background of epicadmium neutrons in the thermal column. The 
increase of the detector activity with the lithium compounds shows that fast neutrons are produced 
by the lithium compounds. 


With each lithium compound two runs were made with the detectors in 
alternate slots. To correct for the different effects of the cadmium-covered 
detectors on the thermal neutron flux at the source the results for the two 
positions were plotted as shown for lithium fluoride in Fig. 4. Here the acti- 
vities have been normalized to the same pile power, corrected to saturation 
and zero decay, and the background from fast pile neutrons subtracted. The 
points for the two runs fall on two curves, whose corresponding ordinates have 
a constant ratio, which is the normalization factor required to bring the curves 
into coincidence with each other. 
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Fic. 4. Detector activity vs. r° for LiF. Note—abscissas shown should be multiplied by 100. 
Run No. 1—WNearest detector is 17.9 cm. from the source. 
Run No. 2—Nearest detector is 28.0 cm. from the source. 
Ratio of ordinates gives normalization factor to bring the curves into coincidence. 
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In Figs. 5 and 6 are plotted the normalized values of Ar? vs. r after sub- 
traction of the background due to epicadmium pile neutrons and the back- 
ground due to photo-neutrons from the deuterium. The fast neutron emissions 
from the samples were obtained from the areas under the curves and the 
indium detector constant previously measured. Table I lists the compounds 
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Fic. 5. Plot of Ar? vs. r for the compounds not containing deuterium. The area under each 
curve is proportional to the neutron emission from the sample. 
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Fic. 6. Plot of Ar vs. r for deuteroxides and the hydroxide. The area under each curve is 
proportional! to the neutron emission from the sample. 


used and the neutron emission observed together with the number of incident 
tritons as calculated from the thermal neutron flux measured at the surface of the 
shells and the amount of lithium in the samples, using 65 XK 10~*4 cm.” per atom 
of natural lithium for the cross section of the Li(m, a)T reaction. In the case 
of lithium hydroxide the calculation of the number of thermal neutrons cap- 
tured per second involves the scattering cross section of hydrogen atoms bound 
in lithium hydroxide molecules. This is not known exactly as it depends on the 
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chemical binding energy of the molecule and its mass. A value of 50 & 107% 
cm.? per atom was used, which is intermediate between those of pure water 
(46 X 107% cm.?) and paraffin (58 X 107%‘ cm.?) (21). The last three columns 
of Table I show the relative yields calculated on the following assumptions: 
(1) The fast neutrons are produced only by (T, 2) reactions on lithium and 
deuterium; i.e., neutron yields from the oxygen, carbon, and fluorine are assumed 
negligible in comparison. 

(2) The thick target yield of neutrons is inversely proportional to the stopping 
power of the compound for 2.6 Mev. tritons. 


TABLE I 
THICK TARGET YIELDS FOR (T, 2) REACTIONS 


| No. of Fast | w ad ' | Calculated relative 

Cittiaiink | tritons | neutrons Neutrons | ae yield 
ave |per second, |per second,| 10 ee vield . — 
| xioe] x10 |" 9 Roasts 


43.9 | 118 | 38 | | 0 | 100* 
| | j 


“a6 | 2. 0 64.6 
| 31.0 | 47. 5.2 | 56.5 50 
| er 1. aa t ae 9.§ 79.8 
| 41.7 | 98.8 | 4 | 91.1 
* Calculated value normalized to 100 for LiF. 
T Calculated value normalized to 46 for LiOD.D20. 


The atomic stopping powers of the light elements relative to air were taken 
from Livingston and Bethe’s review article (11) and the values interpolated 
for the elements not listed. 


Table II shows the neutron yields obtained by Roberts (13) for the bom- 
bardment of the same elements with polonium a-particles of 5.3 Mev. energy. 
As these yields are much less than those listed in Table I, (a, 2) reactions cannot 
account for the fast neutrons observed in the thermal column in the presence 
of lithium compounds. Therefore the fast neutrons are attributed to (T, 2) 
reactions of the type shown in Table III. For the following reasons it is be- 
lieved that the neutrons observed are chiefly from the reactions of tritons on 


TABLE II 
THICK TARGET YIELDS FOR (@, m) REACTIONS (13) 


| Target | 


: N s 7 a- ws 
material Neutrons/10’ a-particles 


Reaction Mev. 


Lis(a, n)B? —3.96 
Lit(a, n)B" —2.78 \ LisCOs } 6.6 
O(a, n)Ne? | 12 O: 0.6 
Ca, n)O" —8.4 Cc 1.0 
Fi9(q. n)Na® CF; 75 

CaF 2 50 
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TABLE III 





Neutron producing : Competing particle 
reactions Mev. reactions (exothermic) 


D(T, n)a +17 


Li®(T, 2)Be® 
bata - +16 

Li®(T, d)Li? 

Li6(T, 2)2a +16 

Li®(T, a)He® Li®(T, p)Li8 
La+n +13.6 I| 


Li?(T, 2)Be® +10.5 i 
Li’(T, an)He® ‘ 
+10.5 


La a 
Li?(T, 2”) Be’ Li’(T, a)He® 


ba+ +8.9 
Li?(T, 2n)2a +8.9 


F19(T, 2)Ne*! +12.1 F(T, a)O8 
F19(T, 2n)Ne” +4.6 F19(T, p)F#! 


O'4(T, a)N 
0'4(T, p)O# 
C#(T, pC 
CAT, a)B" 


O«(T, 2) Fs +1.5 





C(T, n)N™ +4.1 





H1(T, n)He® —0.76* 

















* Threshold = 3.04 Mev. for tritons as bombarding particle (18). 


lithium and tritons on deuterium with contributions of less than 10% of the 
observed yields due to the other elements in the compounds: 


(1) For all the substances used except the deuterides, the observed relative 
yields agree with the relative yields calculated on the assumption that all the 
neutrons are from the reaction of tritons on lithium. 


(2) The neutron yield per triton for the lithium deuteroxide paste is greater 
than for the hydroxide and is much greater than the yield, relative to the other 
compounds, expected from the T + Li reaction. The increased yield is be- 
lieved to be caused by neutrons from the reaction D(T, m)a. 


(3) In Fig. 7 the curves of Ar? vs. r are normalized to the same area, i.e., the 
same number of neutrons from the source. Within experimental error the points 
for lithium carbonate and lithium fluoride fall along a common curve, indi- 
cating that the neutrons from these compounds have the same energy dis- 
tribution and are from the (T, m) reaction on lithium which is the only element 
common to both substances. In the case of lithium hydroxide the large slowing 
down power* of the hydrogen affects the distribution of the indium resonance 
neutrons about the source and alters the shape of the distribution curve. 


* Slowing down power = (scattering cross section per unit volume) X (average logarithmic energy 
loss). 
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Quantitative chawtatitin for these effects are difficult to calculate, and the 
slowing down distribution for the hydroxide cannot be compared exactly to 
the others. Qualitatively the presence of hydrogen is expected to decrease the 
slowing down length of the fast neutrons and therefore shift the distribution 
curve towards the source, as shown in Fig. 7. 














20 





40 60 
DISTANCE FROM SOURCE (CM) 






Fic. 7. Plot of Ar* vs. r normalized to equal area under the curves. The curve for the compound 
containing deuterium lies higher at large distances from the source than the curve for the compounds 
without deuterium, which indicates that the neutrons from the T(d, n)a reaction are more energetic 
than those from the (T, n) reactions on lithium. 








(4) Table III shows that most of the possible (T, 2) reactions on lithium lead 
to three product particles. This means that, depending on the angular distri- 
bution of the emitted particles, the neutrons have any energy between zero 
and the maximum possible. Therefore the mean energy of the neutrons from 
the Li (T, 2) reactions is expected to be less than the maximum energy cal- 
culated from the Q value. This is known to be so for the Li’ (d, m) reactions, 
which show resonances at deuteron energies of 0.65 and 1.02 Mev., corres- 
ponding to excitation levels of 17.17 and 17.45 in the Be® nucleus. The Li® 
(T, m) reactions also form the compound nucleus Be® excited to nearly the 
same energy (17.7 Mev.+ 2/3 of bombarding energy) as calculated from the 
masses and can be expected to produce a similar neutron distribution. There- 
fore the average energy of the neutrons from tritons on lithium is expected to 
be less than that of the neutrons from the T(d, 2)a reaction, which emits only 
two particles and has a higher Q value than any of the lithium reactions shown 
in Table III. This expectation is supported by Fig. 7, which shows that the curve 
for the deuteroxide lies above those for the other compounds at large distances 
from the source, indicating that the deuteroxide yields more energetic neutrons, 
which have a longer slowing down length than the lithium neutrons. 




















(5) Table III shows that, with the exception of hydrogen, (T, m) reactions 
are energetically possible on all the elements present in the compounds used. 
However, it is believed that if the carbon, oxygen, or fluorine contributed more 
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than 10% of the total yield a noticeable difference in the shape of the distri- 
bution curves would appear unless the neutrons produced in each case have 
similar energy spectra. This similarity is unlikely in view of the very different 
Q values for these three elements. Moreover, the relative yields of neutrons 
from LizCO;, LiF, and LiOH.H:O would not agree with the relative yields 
calculated on the assumption that all the neutrons are due to lithium (Table I) 
if there were appreciable contributions to the neutron yield from the other 
elements. 


Average Cross Sections 


The atomic cross section, ¢, for the production of neutrons by bombardment 
of lithium with tritons is a function of the triton energy, i.e., of the residual 
range of the triton. The thick target* yield of neutrons produced by 2.6 Mev. 
tritons may be expressed as follows: 


0 
No. of neutrons per triton = Nj; % a(x)dx 
t 
= N,i¢wR:, 
where V;; = number of atoms of lithium per cubic centimeter of target, 
« = cross section for neutron production per atom of lithium, 
= a function of the residual range of the triton, 
x = residual range of the tritons in the target, 
R; = range of 2.6 Mev. tritons in the target, 
Tay = average atomic cross section for 2.6 Mev. tritons. 
Since the range of the tritons is inversely proportional to the stopping power 
of the target 
_(Stopping power of air) 


we have R; = R, — =< _ = R,(N,S./N: Si), 
(Stopping power of target) 





where R, = range in air at S.T.P. of 2.6 Mev. tritons, 
Sa = molecular stopping power of air, 
S: = molecular stopping power of target, 
N. = Number of molecules per cubic centimeter of air at S.T.P, 
N; = Number of molecules per cubic centimeter of target. 


We have therefore: 
(Na Sa) I 


No. of neutrons per triton = Ny; oa, Ra 
(MN: S:) 


Substituting the known values: 
Range of tritons in air (6)= R,= 6.0 cm. at 15°C. and 760 mm. press., 
Molecular stopping power of air (11)= S,= 2, 
Number of molecules per cubic centimeter of air = Na= 2.55 X 10" at 
15°C. and 760 mm. press., 


* Thick target: target in which the tritons are completely stopped. 
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we obtain: 










i Number of neutrons per triton = 30.6 X 19 Mui) Cline 

it (Ni S2) 

1s N 

ls For a pure chemical —“ is the number of lithium atoms per molecule and is 

[ t 

7 independent of the density or physical form of the compound. For the powders 
used about half of the volume was air but as the stopping power of the solid 
is about 1000 times as great as air the tritons expend most of their energy in the 
solid and a negligible amount in the air. Table IV shows the average cross 

t section for neutron production by 2.6 Mev. tritons as calculated from the 





observed yields. The mean value is 1.5 X 107% cm.? for lithium. 







TABLE IV 


| Nui 


















No. of neutrons i Fan, 

Substance | per 10° tritons nt pole cm.?X 10% 
LiF | 26.8 18.3 1.46 
LizCO; 17.2 | 11.8 1.45 
LiOH.H:0 | 15. | 9.11 1.67 












Mean value Gav= 1.5 X 10-*cm?? 








Using the average cross section given in Table IV we can calculate the yield 
of neutrons from the reaction of tritons on lithium in the deuterides and 
obtain: 
Observed yield for LIOD.D,O + D,O = 21.4 X 107 neutrons per triton 
Estimated yield of tritons on lithium = 9.3 X 107° neutrons per triton 
12.1 X 107 neutrons per triton 
yield due to tritons on deuterium. 
Observed yield for LiOD = 23.7 X 10~* neutrons per triton. 
Estimated yield of tritons on lithium = 20.3 X 107 neutrons per triton 
Difference = 3.4 X 1075 neutrons per triton 

= yield due to tritons on deuterium. 








Difference 






ll 










Using these yields we obtain the average cross section of 2.6 Mev. tritons 
on deuterium: 






LiOD.D.O + D.O o,g,= 0.41 X 107% cm.? 
LiOD Cax= 0.31 X 10-* cm.? 







The value for LiOD depends on the difference of two nearly equal quan- 
tities (observed yield — estimated yield of tritons on lithium) and is much 
less accurate than the value from LiOD.D.O + D.O where this difference is 


larger. 








The D(T, 1)a cross section can be compared to the cross section measured 
by Bretscher and French (5), Baker et a/. (4) and Taschek et al. (17) up to 3.7 
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Mev. triton energy. By using the range-energy relation for tritons (11) and 
0 


their results, ¢ was plotted against x and the integral | o(x)dx = 3.0 X 10-*4 
R; 


cm.” obtained graphically for 2.6 Mev. tritons. This gives: 


—24 
Gay = KO we 0.50 X 10 cm? 


6.0 


which may be compared with 0.41 X 107~*4 cm.* obtained in this experiment. 


It appears that with 2.6 Mev. tritons the average cross section for neutron 
production is greater with lithium than with deuterium. This in part may 
be due to the fact that some of the lithium reactions can produce two 
neutrons per disintegration. For a thick target of pure lithium the yield would 
be one neutron per 1130 incident tritons of 2.6 Mev. energy. 


Slowing Down Lengths 


In Fig. 7 the curves of Ar? vs. r have been normalized to the same area and 
are plotted together with the curve obtained with a Ra-a-Be source in the 
same graphite. 


Assuming a spherical shell source, corrections can be calculated to obtain 
the distribution from a point source (20). The activity at a distance r from 
the shell source is multiplied by a factor (exp, F\/( sinh =) to obtain 

40 rR 20 

the activity at the same distance from a point source with the same total 
emission. R = radius of shell, @ = L?= symbolic age, r = distance from center 
of source. The symbolic age obtained from the uncorrected curves was used to 
calculate the correction factor to be applied to convert the readings to a point 
source. The corrections are less than 2% up to 60 cm., 4% at 80 cm., and 7% 
at 100 cm. Calculations by Mr. S. Kushneriuk also showed that the corrections 
due to lack of symmetry in the source are small. 


Examination of Fig. 7 shows that the curve for (T + Li) neutrons lies lower 
at large values of r than the curve for the Ra-a-Be source, indicating that there 
are fewer very high energy neutrons from (T + Li) than from the Ra-a-Be. 
The curve for LiOD.D.O + D.O lies higher at large ry than any of the others, 
indicating that the neutrons from this are more energetic. The slowing down 
lengths L, obtained are shown in Table V. These were obtained by graphical 
integration of the expression: 


L? -i/ Arr | | Arar. 
0 


0 


About 20% of the integral of the Ar‘ curve is due to points beyond 80 cm. 
from the source where the statistical accuracy was poor. The common curve 
shown in Fig. 7 beyond 80 cm. from the source for the lithium compounds not 
containing deuterium was used to extrapolate the Ar‘ curves to infinity in the 
cases of LiF, LixCO;, and LiOH.H.2O. Subtracting the effect due to T + Li 









and 


24 
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neutrons from the curves for the LiOD.D,O + D,0 gives the slowing down 
length for T(d, a)m neutrons = 23.4 cm. which is much less than expected from 
the elementary age-velocity theory (about 33 cm. for 14 Mev. neutrons). 
Owing to the low intensity at large distances from the source the measured 
value of slowing down length for the T(d, a)m neutrons may be considerably 
in error, but it is believed to be significantly lower than the calculated value. 
This indicates that inelastic scattering of the fast neutrons occurs in the 
graphite so that fewer collisions are necessary to degrade their energy to indium 
resonance than would be expected assuming elastic scattering. 














TABLE V 
| co © 
Ar'd Ard 
5 [ |. sas I. ee Slowing down length (cm.) 
Target Counts/min. in graphite 
(Density 1.64 gm./cc.) 
X 10-" cm.5 | x 10-* cm? 

LiF 33.0 | 14.8 19.3 
Li,CO; 19.6 9.01 19.0 
LiOH.H,O 12.3 5.91 18.6 
LiOD.D,O 30.9 10.96 21.6 
LiOD 28.7 12.4 19.7 
Ra-a-Be 5.73 2.35 20.2 














Errors 


The counting rates of the indium resonance detectors were such that the 
statistical standard deviation was less than 2.5% out to r = 50 cm., equal to 
10% at r = 80 cm., 20% at 90 cm., and about 35% at 100 cm. These are per 
cent errors in the activity due to neutrons from the source after background 
is subtracted. Since only about one-third of the area under the Ar? curves is 
beyond 50 cm. and only one-eighth beyond 80 cm., the error in the area due 
to statistics is about 5%. 

The greatest uncertainty in the experiment is in calculating the number of 
neutrons captured by the lithium from the neutron flux measured at the surface 
of the shells by means of manganese foils. The flux as measured by the foils 
depends on the scattering properties of the shells as well as on the absorption 
cross section and cannot easily be related to the number of neutrons captured 
by the lithium. Approximate calculations, carried out by Mr. Kushneriuk, 
are believed to be accurate to about 10%. 


Owing to the flux gradient in the thermal column the thermal neutron dis- 
tribution is not symmetrical about the source, and hence the shell source is 
not symmetrical. Calculations by Mr. Kushneriuk have shown that the yield 
obtained from the Ar? vs. r curves due to this are 5% low. This has been 
corrected for. 


The absolute yields depend on the value of the neutron emission from the 
Ra-a-Be source which was measured in terms of two standards and is believed 
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to be known to 5%. However the method of comparing this yield with the 
yields of the lithium compounds by the spatial integral of the resonance acti- 
vities is not believed to give better than 10% accuracy, owing to the difficulty 
of locating the foils and source accurately in the thermal column and the poor 
statistics at large distances from the source and the high background of fast 
neutrons in the thermal column. 

The absolute value of the thermal neutron flux as determined by the man- 
ganese flux monitor at the face of the shell is estimated to be accurate to 
about 10%. 

Considering all errors the absolute yields are believed accurate to 20%. 
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ON THE DIFFERENTIAL EQUATIONS OF DIFFUSION! 







By J. D. BABBITT 






Abstract 


The fundamental bases of the differential equation of diffusion are examined. 
From a dynamical equation defining the motion of the gas, an equation of con- 
tinuity expressing the law of conservation of mass, and an equation of state giving 
the relation between concentration and pressure, the differential equations are 
derived for the interdiffusion of two gases, for the diffusion of vapors, and for the 
diffusion of gases and vapors through solids. For the diffusion of gases through 
adsorbing solids, the dynamical equation of the flow is obtained by equating the 
space derivative of the spreading pressure of the adsorbed film to a resistive force 
a to the product of the coefficient of resistance and the velocity of the film. 
The differential equations derived on this assumption agree qualitatively with 
measurements for the diffusion of gases through metals when the adsorption can 
be represented by Langmuir’s equation. When the adsorption follows the BET 
equation, qualitative agreement is found with the diffusion of water vapor 
through hygroscopic materials. It is also shown that Fick’s law is not generally 
valid as the fundamental equation of diffusion. 


















Introduction 







During recent years, the flow of gases and vapors through solid materials 
has been extensively studied. The solids investigated have ranged from metals 
to materials of high porosity, and an attempt has been made to explain the 
flow by a great variety of mechanisms. Where the flow has been through finite 
pores within the solid and where there has been an absolute pressure difference 
so that the flow has been hydrodynamical rather than diffusive, the amount 
of gas transferred has been satisfactorily explained over all ranges of pressures 
and degrees of porosity by various theoretical formulas. Where, however, 
the flow has been diffusive, either the interdiffusion of two gases through a 
porous barrier or the diffusion of a single gas into a nonporous solid such as 
a’metal, the results have not been satisfactorily explained, and, in consequence, 
it is desirable to re-examine the fundamental assumptions of diffusion. 












Whenever a theory of the diffusion of gases through solids is presented, 
Fick’s law is invariably used as the fundamental equation. Barrer (3) gives 
this law in two forms, first, as the equation for the steady state of flow, 






(1) 





0 
os ees 
Ox 


P 





where the rate of permeation P is given in terms of the concentration gradient, 






Oc : ; ee ‘ 
— , and the diffusion coefficient, D., and, second, in the form, 
Ox 





0°c 
Ox? 


0 
= & 
ot 


(2) 
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which gives the variation of the concentration with time and which is analogous 
to Fourier’s equation for the flow of heat. 


These equations have long been used with great success to describe the 
interdiffusion of gases and solutions. With both gases and solutions, however, 
an assumption can safely be made which does not necessarily hold when a solid 
is present, namely, that the concentration is linearly proportional to the pres- 
sure. When this assumption is true, the substitution of p for c in Equation (1) 
results only in a change of the numerical value of the constant D, and it is not 
necessary to distinguish carefully between the part played by the pressure and 
the concentration. When, however, the diffusion is through solids the relation 
between pressure and concentration is seldom linear owing to the complexities 
of adsorption, and it is essential to know whether Equation (1) or the equi- 
valent equation in terms of the pressure gradient is the correct statement 
of the law governing the diffusion. 


It is very probable that the tendency to use Equations (1) and (2) as the 
fundamental equations for diffusion without examining explicitly the assump- 
tions involved arises from a confusion between the kinetic picture of diffusion 
represented by Meyer’s equation and the macroscopic picture represented 
by a differential equation such as (2). Although in Meyer’s equation (cf. Jeans: 
The Dynamical Theory of Gases—Chapter XIII) the number of molecules 
transferred is proportional to the concentration gradient, it does not auto- 
matically follow that, in an equation such as (2), which is essentially a field 
equation, the concentration is the correct quantity to use. A study of the 
analogous differential equations for the flow of heat, electricity, and magnetism 
shows that the fundamental equation of flow expresses the flux as proportional 
to the gradient of a potential function (the temperature, electromotive force, 
and the magnetic potential, respectively). In diffusion, the corresponding po- 
tential function is the pressure and not the concentration and from analogy 
one should infer that the fundamental equation should be expressed in terms 





of the pressure. 


A second source of confusion has been the failure to distinguish diffusion 
from other flow mechanisms. The term diffusion has been used so loosely that 
it is often difficult to know what is precisely meant. In this paper, we shall 
restrict the term diffusion to two mechanisms of flow corresponding to two 
distinct kinds of resistive forces. Where the flow is confined to the gaseous 
phase, we shall use the term diffusion only when the total pressure is every- 
where equal and only differences of partial pressure exist. Under these con- 
ditions, we assume that the force resisting the movement of gas A through gas B 
is due predominantly to the interference which the molecules of gas B offer to 
the movement of the molecules of gas A, and conversely that the force resisting 
the movement of gas B is due predominantly to the interference of the mole- 
cules of A. When, on the other hand, the flow is through an adsorbing solid, 
gas may be transported in the adsorbed phase, and we shall also apply the 
term diffusion to such transport when there is a movement of the gas in the 
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adsorbed phase from points of high to those of low concentration. Restricting 
ourselves to these two specific mechanisms of flow we shall proceed to derive 
the differential equation which shall govern such diffusion. 


In deriving the differential equation, we are indebted to Muskat (9) who 
was faced with a similar problem in studying the flow of gases and liquids 
through sands. There are differences between Muskat’s problem and ours. 
We are limiting ourselves to diffusion in the restricted sense outlined above. 
Muskat was interested only in flow with an absolute pressure gradient and he 
assumed that his flow could be defined by Darcy’s equation. When these dif- 
ferences are kept in mind, however, Muskat’s problem is exactly the same as 
ours, namely, to derive the partial differential equation for the flow, and it is 
possible to adapt his analysis to our problem. Therefore, to show what is 
involved, we shall proceed to outline very briefly the treatment used by Muskat 
and shall then use this method as a model to derive the differential equations 
of diffusion. 


Muskat’s Development of the Equations for the Flow of a 
Homogeneous Fluid Through Porous Media 


In developing the equations for the flow of a homogeneous fluid through 
porous media, Muskat has pointed out that three separate conditions must be 
satisfied. These conditions are expressed by: 

(i) The dynamical equation of motion, 
(ii) The equation of continuity, 
and (iii) The equation of state. 


The dynamical equation defines the nature of the flow, the equation of con- 
tinuity is the mathematical expression of the law of conservation of mass, and 
the equation of state expresses the relation between the density and the pres- 
sure of the gas. 


To develop the dynamical equation we know that, in general, a volume 
element of fluid will be subject to three types of forces. These are (1) the pres- 
sure gradients of components z, op and - (2) the external ‘‘body 

ox oy 02 
forces’, such as gravity, of components F,, F,, and F, acting on unit volume 
of the fluid, and (3) the forces opposing the motion of the fluid, which are due 
to the internal resistance of friction and which may be written in the form 
A,, Ay, and A,. By equating these forces to the product of the mass and 
acceleration of the volume element, we obtain the generalized equation of 
motion. Thus, in one dimension, 

Du 


op 
cane, Perit eee Seen F, Ax 
. Dt Ox * * 3) 


‘ D : Meee ; 
where, following common usage, D denotes a differentiation following the 
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motion of the fluid, p is the density and, u the component of velocity in the 
x-direction. (Note: Throughout this paper we shall, for convenience, write all 
formulas in one dimension only. The expansion to three dimensions will be 
readily apparent.) 

As Muskat points out, when the velocity of flow is small, the inertial or 
accelerative forces will be negligible compared to the other forces, so that the 
term on the left of (3) may be taken as zero. Although the effect of the gravi- 
tational forces will be appreciable in the interdiffusion of two gases of different 
densities, we are not, at present, interested in that aspect of the problem and 
we shall neglect all external forces. With these assumptions the dynamical 
equation reduces to 


#_ A, = 0, (4) 
Ox 


which states simply that the pressure gradient is balanced by the internal 
resistance or frictional forces. 


For the type of flow in which he was interested, Muskat used Darcy’s law 
as his fundamental dynamical equation. Darcy’s law states that the velocity 
of a fluid flowing through a porous medium is directly proportional to the 
pressure gradient acting on the fluid. Mathematically this can be expressed 


0 

oe 

Ox k 
where uy is the viscosity of the fluid and k is a factor depending on the per- 
meability of the medium. Thus in Muskat’s treatment 


Agim 


k 


The equation of continuity that applies to flow through porous bodies is 
that commonly used in hydrodynamics, namely, 


‘ 0 0 
div (ps) = © (pv.) + © (o0,) + - (pv,) = (7) 


except that it is necessary to introduce the porosity f to compensate for the 
fact that the area of flow has been taken as the area of medium plus pores, 
whereas the volume containing the gas is that of the pores alone. Introducing 
the porosity and considering flow in one dimension only, we get 


0 Op 
ae at) Laces 8 
= (pu) f (8) 


Muskat takes as his equation of state the general relation 
p = po pme??, 


where po is the density at unit pressure and m and Bare constants. With this 
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equation the particular fluids of physical significance may be classified as 
follows: 
Liquids: 
incompressible liquids, 8 = 0, 
compressible liquids, 6 ¥ 0, 
Gases: B=0 


isothermal expansion m = 1, 
: ; : Spec. heat at const. vol. 

adiabatic expansion m = ————_______ 

Spec. heat at const. pres. 

In the work that follows we shall be interested in gases only so that throughout 


we shall take 8 = 0. 


By combining the dynamical equation with the equations of continuity and 
of state, Muskat obtains for the differential equation of flow, the expression 
ett 
Op k Op ™ 
ee ee | (9) 
ot (1 + m) fupo ’” ox 


This equation occupies in Muskat’s theory of the flow of fluids the same posi- 
tion that is occupied in the theory of the flow of heat by the equation 


a6 «8 (10) 


at ax?’ 
where @ is the temperature and « the thermal diffusivity. Muskat goes on to 
point out the implications of this equation and, in particular, to derive the 
expression for the pressure gradient in steady flow and to show the difference 
between isothermal and adiabatic flow. We are not, however, directly inter- 
ested in flow obeying Darcy’s law, and having derived this equation as a model 
for diffusive problems we shall proceed no further with it. 


The Differential Equation for the Interdiffusion of Two Gases 


In applying this analysis to diffusion, we find that the fundamental dynam- 
ical equation for diffusion was developed by Stefan and Maxwell in the nine- 
teenth century but that, in the meantime, the hydrodynamical approach to 
diffusion represented by such an equation has been neglected. The best pre- 
sentation is contained in the article on ‘‘Diffusion” by G. H. Bryan and A. W. 
Porter in the fourteenth edition of the Encyclopedia Britannica and we are 
largely indebted to them for the ideas used in this section. 


The fundamental concept of this picture is that, in the interdiffusion of two 
gases A and B, the portion of the gas A contained in a small element of volume 
v will experience from the gas B a resistance to its motion equal to Cpip2v(u1 — 
us), where pi, p2, and 14, us are the densities and velocities of the gases A and B 
respectively and C is a quantity called the coefficient of resistance. It is ap- 
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parent that the gas B contained in the same volume element will experience 
an equal and opposite resistance Cpip2v(u2— 11). 


With this concept the dynamical equations can easily be set up. Thus 
Maxwell gives for the gas A the equation 


mee + ts Xpi + Cpips (ui — U2) = O, (11) 
Dt Ox 
where /; is the partial pressure of A and Xp; the external forces. It is easy to 
show as before that the effects of inertia and of external forces are negligible 
so that we get, 
Op; 


ax 


+ Cpip2 (ui — U2) = 0. 


(12) 


Similarly for B ue + Cpip2 (u2 — U1) = 0. | 
x 


These are the fundamental dynamical equations for the interdiffusion of two 
gases. 


The total pressure of the system, P = 1+ p2, must remain constant and 
it follows that the fluxes of partial pressure due to the two gases must every- 
where be equal and opposite so that piu1+ pou2= 0. Also, since in diffusion 
all the motions are very slow compared with those set up by mechanical actions, 
we can assume that the flow is isothermal and can take as our equation of 
state that p1= kip: where ki is a constant. Thus 


Ripits + Repetto = 0. 


. U1 U2 Ui — the Ui — ue 
Now —_=-— = 


ps ‘i ek P 


= —e SE ’ 


so that 


Pipe(ui — U2) _ Rikepip2(ui — U2) 
“> => Un = = e 
Pit Pe P P 





Substituting in the dynamical equations, we get 


dpi , CP dp. , CP 
pe Ee = 0 d —— qk2) = 0. 13 
ra i (pitts) an —e. (Pott) (13) 


Whence, substituting 1 = kip: and p2 = kop: and transposing we obtain 


_ Rik Op2_ 


14 
CP dx am 


poke = 


If we put equal to D, the coefficient of diffusion, we obtain Fick’s law 


as expressed in Equation (1). We see, therefore, that (1) is a legitimate expres- 
sion for the interdiffusion of two gases. 
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Substituting Equations (14) into the equation of continuity 


Apu) _ _ 9 
ox at’ 


we get the differential equation 


(15) 


which is equivalent to Equation (2). 

It is apparent, therefore, that according to this picture Fick's law is a cor- 
rect and adequate representation of the interdiffusion of two gases. It should 
not be overlooked, however, that even here the fundamental dynamical equ- 


ation, 
Opi 


+ Cpip2(ui — uw) = 0 
dx 


is expressed in terms of the pressure gradient and it has been possible to trans- 
form this into the commonly used concentration form of the law only by 
the assumption of Boyle’s law as the equation of state. 


The Diffusion of a Condensing Vapor 


It has been pointed out by Stefan (13) that, where the diffusion is that of 
a condensing or evaporating vapor, there is no transfer of the second gas. 
Consider the evaporation of a gas from a liquid surface in a cylinder. Let p; 
and p2 be the partial pressures and let the suffix 1 apply to the vapor and 2 
to the permanent gas. Since there is no movement of the permanent gas 
u2= 0. As before we assume Boyle's law as the equation of state and put 


pi= kip, and po= kop. 


The fundamental dynamical equations are 
Opi 


~ Apo 
-+ Cpipou, =0 and Ops 
Ox Ox 


— Cpipot, = 0. 


Since pi + po = P, these reduce to 


om + ee (pit1) = 0 and > Ae 


ui) = 0. 17) 
Ox kiks Ox kike (Pr ( 


Introducing the equation of state, we get 


en kike P Opi ; (18) 
C(P — p1) Ox 


pill, = 


If we follow Stefan and put aus = D,, the coefficient of diffusion, we have 


pill 
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This is the equation for the mass transferred corresponding to Fick’s law in 
the form (1). Introducing the equation of continuity, we get 
Op, D, 0" 
— = —— — [log (P — p,)). (20) 
ot ky Ox 
Thus we see that, in the diffusion of a vapor, the function log (P — kip:) 
replaces the density p; used in the interdiffusion of two gases. 


It is instructive to study how this equation works in practice. To determine 
the rate of diffusion of a vapor the common procedure is to measure the mass 
transferred from a region of high partial pressure, which we shall denote by 
ip1, to one of low pressure, which we shall denote by 2f;. The equation for 
the mass transfer may be written 


D, @ 
pil, = — — [log (P — p,)]. 
ky Ox 
Since by the equation of continuity piu; is constant when the flow is steady, 
we can integrate this equation. Hence 
D, | log (P= 11) | (21) 


puy= — — — 


hE ot mae 


where L is the distance between the points where pi1= ip: and p; = 2p. We 
can write this equation 


(141 = — De if log ¢ a iPr Be ) ’ (22) 
bi LG P — spi 
and expanding the logarithmic term 
De —» 1 —opi\ ‘ 
pai fbma (mm, ay 
ki LLP—opi 2 \P— opi 


In most experiments on vapors the partial pressures p; will be much smaller 
than the total pressure P and all terms in the series except the first may be 
neglected. Also, in most diffusion experiments 2/; is zero or very close to it so 
that 2p; may be neglected in comparison with P. Our equation thus becomes 


D y 1 are ‘ 4 
Pilly = — Pr ~ 2P 1, (24) 
ki P L 
or, if we express it in terms of the densities 
D, — >» 
pil, = — 6 tea (25) 
P L 


If we write D = $ this equation is identical with the integrated form of the 


equation for the interdiffusion of two gases and we see that, when the partial 
pressure of the vapor is small in comparison with the total pressure, the error 
introduced by using Fick’s law in the form 
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Rik» Opi 


pm) = — —— 


CP dx 


Op. 


Ox 


is negligible. 


Diffusion Through an Inert Porous Medium 


It is obvious that when we have two gases separated by a porous barrier the 
resistance to the mixing of the gases will be greater than if the barrier were not 
present. In general, there will be two effects: 


(i) Since the area available for diffusion will be much less and since, in general, 
owing to the tortuosity of the paths, the molecules will have to diffuse a longer 
distance, the rate of diffusion will be less than if the barrier were not present. 
To represent this reduction, we can assume that the diffusion coefficient is 
reduced proportionally to the reduction in area (which, in turn, is proportional 
to the porosity f of the medium) and also that it is reduced by a factor r = 


l : ; Pe 
— = ratio of length of barrier to effective path traveled by the molecules. 


e 


Hence the equation of mass transfer becomes 


pil, = ~~ ° (26) 


(ii) In developing the equations for diffusion it has been assumed that the 
resistance to the movement of the gases is the result of the mutual interference 
of the two gases only, and that, in particular, this resistance is so great and 
the movement of the gases is so slow that any hydrodynamic or viscous resis- 
tance can be neglected. When we are considering the diffusion through porous 
bodies it is not immediately obvious that these assumptions will hold, for if 
the pores are very small it is conceivable that the resistance to hydrodynamic 
flow might be comparable to the diffusive resistance. According to a kinetic 
picture, this would take place when the collisions of the molecules with the 
walls of the pores became comparable in number to the collisions with mole- 
cules of the second gas. For the present, however, we shall neglect any hy- 
drodynamic or viscous resistance and assume that all resistance to the flow 
is diffusive. The measurements of Penman (10) indicate that this assumption 
is valid, at least for diffusion through soils. The theoretical justification for 
the assumption is that at atmospheric pressure the size of the pores would 
have to be very small (< 1075 cm.) to be comparable to the mean free path of 


the molecules. 
Writing Equation (26) in the form 


0 
pil, =_l po 
Ox 
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we take it as our equation of transfer for the diffusion through an inert porous 
medium and we can proceed, as for simple diffusion, to develop the differential 
equation. 
The equation of continuity is similar to that used by Muskat for the hydro- 
dynamic flow through such media, namely, 
i 
Ox Ot 


Combining this with Equation (27) we obtain 


Opr _ D’ d°p1 
ot Ox? 

or, 9p _ rD O"p1 | (28) 
ot Ox? 


Thus the differential equation is similar to that for the interdiffusion of two 
gases except for the constant factor r. This result is confirmed by the measure- 
ments of Penman (10). 


Diffusion of a Vapor Through an Inert Porous Medium 


It is obvious that, in the range of vapor pressures where the approximation 
developed above is applicable, the diffusion of a vapor through an inert porous 
medium will follow the same laws as the interdiffusion of two gases. 


Diffusion of a Vapor through an Adsorbing Medium that is not 
Permeable to Nonadsorbed Gases 


In discussing the movement of a gas through a medium by which it is 
adsorbed, we shall restrict ourselves, in this section, to consideration of a non- 
porous medium in the sense that there are no channels or pores of sufficient 
size to allow the penetration of a nonadsorbed gas. It is necessary, therefore, 
to conceive of the gas as being adsorbed on the external surface and migrating 
through the medium as an adsorbed film on internal surfaces. There is con- 
siderable supporting evidence for such a picture and it is not necessary to 
make any assumptions about the internal surfaces save that they are accessible 
to adsorbed molecules and nonaccessible to the gas phase. We shall be limited 
only by the assumption that, once adsorbed, the gas migrates as an adsorbed 
atom or molecule until it reaches the external surface on the low pressure side 
of the barrier where it evaporates. 


It is obvious that, on such a picture, some of the fundamental ideas of dif- 
fusion must be modified. In the interdiffusion of two gases, the basic assumption 
is that the resistance to the movement of the two gases is due to the relative 
motions of the two kinds of molecules, and that this relative motion gives rise 
to a resistive force proportional to the densities of the two gases and to their 
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relative velocities. In the diffusion of an adsorbed gas, as pictured above, no 
second gas is present within the adsorbent to resist the motion of the diffusing 
gas, and it is only possible to conceive of the resistance as residing in the 
“friction” between solid and adsorbed film. In the modern treatment of ad- 
sorption, such a resistance is related to the energy required to move an atom 
over the potential barriers between adsorbing sites on the surface. 


It is evident, according to this picture, that the second gas, if present, plays 
no function in the diffusion save that of hindering the movement of the dif- 
fusing gas up to the external surface of the medium. Since the resistance thus 
offered is in general, negligible compared to the resistance within the medium, 
there will be no distinction here between hydrodynamical flow, where there is 
only one gas, and diffusion, where a permanent nondiffusing gas is present in 
addition to the adsorbed gas. In such a diffusion, it is immaterial whether the 
pressure difference of the diffusing gas is a total or a partial pressure difference. 


To find the differential equation governing the movement of a gas under 
these considerations we return to our fundamental dynamical equation which 
we express in the simplified form 


sen 
Ox 


The problem is (1) what is the pressure gradient that applies to such a diffusion 
and (2) how are the resistive forces A, related to the physical characteristics 
of the system. 


In the statistical treatment of adsorption, an adsorbed layer is found to have 
the properties of a two-dimensional gas. In particular, it is possible to derive 
a potential function ¢ which has, for an adsorbed film, exactly the same thermo- 
dynamic properties that the pressure p has for a three-dimensional gas. ¢ has 
been called the spreading pressure and is related to the free energy of the 


adsorbed film by the relation ¢ = — , where A is the area of the film. 


@ is a measure of the force that must be applied to keep the film from spreading. 
It would seem, therefore, from these qualitative considerations that the spread- 
ing pressure ¢ is the appropriate potential function for the fundamental dynam- 


ical equation. 


In considering the hydrodynamical picture of the movement of an adsorbed 
film, it is necessary to think of the film as moving as a whole over the surface 
with a velocity of motion u just as, in viscous flow, a gas is pictured as moving 
as a whole even though it is known that the viscous forces are due to random 
motions of the molecules. We shall, therefore, treat the adsorbed film as having 
a velocity of motion u relative to the surface and shall assume, as is natural, 
that the resistive force to the motion will be proportional to u. We shall put 
the resistive force equal to Cu where C is again the coefficient of resistance. 
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There is, at first sight, no reason for thinking that the resistive force will be 
independent of the number of molecules adsorbed. To the contrary, in the 
kinetic picture of adsorption, it is sometimes assumed that the molecules 
adsorbed at low pressures are bound more tightly, and consequently would offer 
greater resistance to flow, than the molecules that are adsorbed at higher pres- 
sures. It would seem, therefore, that, ‘rom kinetics, the resistive force might be 
taken as dependent on the number of raolecules adsorbed. On the other hand, 
from the statistical point of view, vari-tion of binding force with concentration 
is already taken care of in the expression for the spreading pressure—when the 
atoms are bound firmly they contribute a small effect to the spreading pressure 
and vice versa—and it is not necessary to introduce the influence of binding 
force into the resistive force. We shall, therefore, express our resistive force as 
Cu and assume, as justified by results, that C is a constant, characteristic of 
the gas and the medium. 


The fundamental dynamical equation is, therefore, 
0 
Sie (29) 
Ox 


where the sign is positive since uw is in the direction of decreasing pressure. 


The equation of state for such a diffusion is given by the relation between @ 
and the concentration of adsorbed molecules. This relation, of course, depends 
on the adsorption isotherm and is readily available only for those isotherms 
that have been statistically derived. We shall not then, at this time, attempt 
to give the generalized differential equation for flow through adsorbing media 
but shall confine ourselves to those types of adsorption for which a statistical 
expression for ¢ is available. 


The Differential Equation for the Diffusion of a Mobile Monolayer 


The simplest model of an adsorbed film on a solid is that designated by 
Fowler and Guggenheim (5) as a mobile monolayer. In a mobile monolayer 
some point in the adsorbed molecule (for example its center of mass) is bound 
tightly to the surface in the normal direction but can move with complete 
freedom in the two directions on the surface. It is unlikely that such an ideal- 
ized picture will represent the behavior of any film actually adsorbed on solids 
but since it is the simplest representation of adsorption we commence with a 


discussion of it. 


For a mobile monolayer, Fowler and Guggenheim give for the spreading 
pressure ¢, the relation 
oe NkT (30) 
A 


where N is the number of atoms adsorbed on a surface of area A. For con- 
venience we take A as the area of surface per unit volume so that N equals 
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number of atoms adsorbed per unit volume. It is apparent from the similarity 
with the equation 
NkT 
ps = 
I 
in three dimensions that a mobile monolayer is the two-dimensional analogue 
of a perfect gas. 


When we substitute this value for ¢ in Equation (29), we get 


~ aT ON (31) 
CA ox 


and for the mass transfer 
T ,,a1 
Nan ~ Se ne. (32) 
CA ox 
Since p, the density of adsorbed gas, is clearly proportional to N, we may 
write the equation of continuity in the form 
0(Nu) es: (33) 


Ox 


which is more convenient when we are considering adsorbed gases. 


From (32) and (33) we obtain the differential equation for the flow, namely, 


aN _ kT 8 (2%) 


‘at CA ax 
= (34) 


Ox 


It is perhaps worth while to look at the implications of this equation. For 
dN 
= 0, we have 


steady-state flow, since 
at 


2 AT2 
en = 6 (35) 
Ox? 


To find how the concentration would vary with distance we consider a 
barrier of thickness L, with concentration No at surface x = 0 and WN; at 
x = L. Integrating (35) and putting in the boundary conditions, we find 


— N?) - + N?. (36) 


Thus the concentration varies as the square root of the distance. 
We find how the diffusion coefficient of the commonly used equation 


ta ee 
Ox 
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might be expected to vary with concentration we write the equation for the 
mass transfer in the form 
kT dN 
“= ——— (38) 
2CA 0x 
and integrate from x = 0 to x = L, remembering that, by the equation of 
continuity, Nu is independent of x. Thus 


2 2 
_ #f Ne Nt, (39) 
2CA L 
Similarly we integrate (37) bearing in mind that P is equivalent to Nu and 
Cto N. Thus 


7. 


Nu = P = D,~°— A. (40) 


‘ ; sts. NE 
Equating and writing —— = D, we get 
CA 


Not Nr. (41) 


D, = D 
2 


Thus we see that, according to this picture of a mobile monolayer, we should 
expect the diffusion coefficient as commonly derived from Fick’s law to be 
proportional to the mean concentration. 


The Differential Equation for the Diffusion of an Ideal Localized 
Monolayer (Langmuir Adsorption) 


According to the model of a localized monolayer developed by Langmuir, 
the molecules are attached to the surface at definite points of attachment 
which are generally referred to as sites. From a kinetic argument, Langmuir 
was able to show that the isotherm for such an adsorption is represented by 
an equation of the form 

ce (42) 
botp 1+ap 


‘ ; ; ; RY 
where @ is the fraction of sites occupied and pp) = — is a constant. 
a 


Fowler has shown how this equation may be statistically derived and has 
obtained an expression for the spreading pressure (5). He gives 


N, 

aie log a's (43) 
A 1-0 

where NV, may be taken as the number of sites, and A the surface area, per unit 

volume. 


FE 


Substituting this value of @ into our fundamental dynamical Equation (29), 
we get 
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NkT 0 ('oe : ) + Cu =0. (44) 
A oe 1—@ 


‘ N : 
Putting 6 = —, where N = number of adsorbed atoms and transposing we 


N N 
ht = ae. Se 45 
F AT 9 (tog —*s) ~ 


find 


Whence the equation of mass transfer can be written 


Nu = — NokT v2 (og ) 
CA Ox N,-N. 


_NiakT NaN. 
CA N.—wN Ox 


On comparing this equation with Fick’s law 
we see that 


(47) 


N,kT ‘ ; 
-—— has been plotted against @, and the resulting curve shows 


/ 
In Fig. 1, D. | 


how the diffusion coefficient of Fick’s law would vary with @. 


Equation (47) and Fig. 1 give the instantaneous value of D, at any concen- 
tration @ but in practice it is always an integrated value that is measured. 
To see how the measured values will vary with concentration we proceed with 


our analysis. 
The equation of continuity is again 
O(Nu)_ _ aN 
ox ot 
and introducing the fundamental dynamical Equation (45), we obtain 
aN _N,kT 9 (5: N on) 
at CA dx \N,—N Ox 
_N.kT 2 
CA ax? 


This is the differential equation for diffusion in an adsorbing solid when the 
adsorption obeys Langmuir’s isotherm. 


[V, -N —N, log (N, — N)). (48) 
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In the steady state, we should have 


S; iN, —N — Ny lee (8) ~ oO, (49) 
7 


and the equation giving the variation of the number of adsorbed atoms with 
the distance x is 


N,—WN,z x 
N,—N-—VN, log (N, -—N =(w —~ N;, —N,1 Me Ne) 2 
og ( ) 0 L ST. N, L 


+ N,— No— Nz log (Ns — No). (50) 


Fic. 1. Variation of the coefficient of the differential equation when the adsorption obeys 
Langmuir’s equation. 


The rate of diffusion of a gas through a solid is most commonly determined 
by measuring the mass transferred through a definite thickness of the solid 
with a constant pressure difference. We shall now calculate the amount that, 
on our picture, would diffuse in such an experiment. 


The equation for the mass transfer can be expressed 


NskT = [V. — N—N, log (N; — N)], 


CA ox 


Nu=-— 


and integrating from x = 0 to x = L, we get 


Nu = = - E — N—N, log (N, — x) | 


7 — 
“ - NAT 1 yy p+, log N, =, 
cA N,—N, 


N=No 


N=Nyz 
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’ 


or if we wish the expression in terms of 6 = 





8 


Mar | =f 
Ne « ~~ 6h . 53 
a kt ee 












This expression gives the amount of gas diffusing through a barrier when 
6 = @ on one surface and 6 = 6; on the other. To find, in practice, a diffusion 
following these laws we turn to the diffusion of gases in metals. It is well known 
that certain gases are adsorbed by metals in such a manner that the isotherm 
can be represented by the Langmuir equation and, according to our picture, 
the diffusion of such gases through the metals should follow the equations 
developed above. 









To test these equations the experimental results of Smithells and Ransley 
(12) are most convenient. Their measurements give the volume of gas at 
N.T.P. diffusing per second through unit area of metal, 1 mm. thick when 
there is zero pressure on one side and pressure p on the other. Since volumes 
at N.T.P. will be proportional to the mass, Smithells and Ransley have essen- 
tially measured Nu as given by Equation (53) and their results may be used 
as a direct test for that equation. 
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Fic. 2. Diffusion of hydrogen through copper. @ Experimental values from Smithells and 
Ransley. 









In Fig. 2 we present the results of Smithells and Ransley for the diffusion 
of hydrogen through copper at 723°K. and in Fig. 3 the results for the diffusion 
- of hydrogen through iron (etched) at 686°K. The results have been expressed 
as a fraction of the weight diffusing at the highest measured pressure and the 
equation is made to agree with experiment at that pressure. As a second 
adjustable constant we have a in the expression 0 = — and the curve 
a 
for the value of a giving the closest fit to the experimental results is shown in 
the figures. At present it is impossible to deprive supporting evidence for the 
theory by comparing the value of a from diffusion with that from the isotherm 
since the isotherms have been determined only at temperatures considerably 
lower than those at which diffusion experiments have been made. It can only 
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be said about a that the value used in Fig. 2 is the same as that considered 
satisfactory by Smithells and Ransley when derived from another theory. 


The agreement with experimental results shown in Figs. 2 and 3 is encour- 
aging but, since it is not the intention in this paper to study in detail the dif- 
fusion through metals and its relation to the equations above, we shall proceed 
no further with this analysis. We wish to point out only that the simple 
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Fic. 3. Diffusion of hydrogen through iron (etched). @ Experimental values from Smithells 
and Ransley. 


assumptions made in deriving our differential equation are sufficient to lead 
to an equation that is, at least, in qualitative agreement with measurements, 
and that the theoretical basis of this equation is more substantial than that 
of the semi-empirical equation that has hitherto been used. 


The Differential Equation for the Diffusion of a Gas Obeying 
the BET Isotherm 


There is one other class of adsorption isotherm that has been represented by 
a statistical equation and that has great importance because of the variety of 
solids that adsorb vapor in that way. This is the equation of Brunauer, Emmett, 
and Teller for polymolecular adsorption. The BET equation was originally 
derived, like Langmuir’s equation, from kinetics but recently Cassie (4) and 
Hill (6) have shown how it may be obtained from statistics, and from their 
analysis it is possible to derive an expression for the spreading pressure ¢. 
We shall, therefore, proceed to discuss how a gas adsorbed by the BET mech- 
anism would, on our theory, diffuse through an adsorbing solid. 


To obtain the expression for the spreading pressure ¢, we have that 


OF ads 


a 


and from Hill we find that 
Fias = — kT log Q, (54) 


where Q is the complete partition function for BET adsorption. Evaluating 
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Q in terms of the total number of adsorbed molecules NV, the number of localized 
sites B, the number of molecules X adsorbed on the localized sites, and the par- 





€1 *, 





tition functions j,e*7 and jze*? for the molecules adsorbed in the first layer 
and in higher layers respectively, we find for F,a, the expression 


Fias = — er | log N + Blog B— (N—X) log (N—X) — (B—X) log (B—-X) 






—2XlogX +X =~. + (N—X) +X log js + (N —X) lozit | (55) 






Whence ? eal ae OF ads ne B OF ads 
dA A 0B 
BkT B 
ane WE mere * 56 






Substituting this value of ¢ in our fundamental dynamical equation (29), 






we get 
BkT 2 ( so ee ) 
—~ = (log ——_ ) + Cum, 
A 0x © BX 
or / “= - FET 2 (e575): (57) 
CA ax B-X 





and the equation for mass transfer becomes 
Nu = _ BkT n2 (tog — 
CA ax B-X 
Re ee (58) 
B-X 0x 







To eliminate X, we use the equation given by Cassie that 
(N — X) (B — X) = BX? 
€r om. Os 


where B= Jue Fri tin (59) 
Js 


































es aX_ ss B-X aN (60) 
ax N+B—2(i —B)X dx 
oe ya N+B — oes 4(1 — 8)NB- 61) 






Following Cassie, we use the minus sign for X as giving the minimum free 
energy, and substituting, we get 
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It is very difficult to observe N experimentally and in practice the measured 


p 


quantity is almost always the pressure p. We shall use the notation P = >. ’ 
0 


where fp is the saturation vapor pressure. In this notation the BET equation 
can be written 
ea = : (63) 
B (1 —P) (1 -—P+ CP) 





where C = 1/8. 
With the BET equation we can eliminate N from (62) and we get 
_ BRT 1 
CA y/ x (48—2)(1—P)(1— P+CP) , (a —P)(1—P+CP)} 
CP Cr 
 C— 4 CP) a 
{1 —P)(l-—-P+CP)}? dx 





Nu = 





(64) 


In experimental work on adsorbing solids it has been usual to express the 
results in terms of an equation 
oP OE 
Mass transfer = Nu = — D,—°* (65) 
Ox 
By a comparison of (65) with (64) we see that D, would be expected to vary 
with P according to the relation 











ee wee oe aS ae con 
"CA /, 4 G-D0-P)G-P+CP) , (d—P)—- P+ CPP 
Vit : += = 
CP cr 
C(1 — P? + CP*) “ 
wa een errr (66) 
(i — Fl ~ Faery 
3 | kTB? : 
In Fig. 4 D, | CA has been plotted against P and the curve shows the 


way the diffusion coefficient of the simple equation (65) will vary with pressure 
if the adsorption follows the BET isotherm and our assumptions are correct. 
In calculating D, from (67) 8 and C have been given the values that would be 


appropriate to spruce wood, that is, C = ; = 11.80 (1). Although spruce wood, 


being permeable to nonadsorbed gases, does not come within the class of 
solids we are considering we have taken this value for C as being readily 
available and as giving a reasonable figure. We are interested at this stage in 
only a qualitative comparison and wish only to see the general behavior. 


In practice it is difficult to measure the differential value of the diffusion 
coefficient as given by (66); nearly always Equation (64) is integrated 
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before it is applied. We shall, therefore, proceed to derive such an integrated 
equation and show how the commonly measured coefficient of diffusion is 
related to the mean pressure. The first step is to obtain the differential equation. 


120 r 


100 + 





2 oh specie 


Q2 O04 06 Os 
RELATIVE HUMIDITY 


Fic. 4. Variation of the coefficient of the differential equation when the adsorption obeys the 
BET equation. 


We may write the equation of mass transfer (62) in the form 


ON 
r (Ox 


(67) 
where Y = (N — B)? + 48NB. (68) 


Equation (67) is equivalent to 


.vre B (28 — 1) log] VY + N + B(26 — »}. (69) 


CA dx 

and using the equation of continuity in the form 
a(Nu) _ _ ON 

Ox at’ 


we obtain for the differential equation 


oF ae mV - B (28 —1) log] VY + N+B(26- v}. (70) 


To find the expression for the integrated mass transfer, we consider the 
movement of a gas through a barrier of adsorbing material from the surface 
x = 0 where P = Po to the surface x = L where P = P;. Integrating Equa- 
tion (69) we get 
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Nu = Bet VY B (28 —1)log | V¥+N+8 (28 — » |} 


N=No 


N=N, 


Similarly integrating Equation (65) we have 


Nae Bones. 


Hence equating (71) and (72) we get for D, 


SVy_ Se y 2s 
_ Ber} B(28 1) log] V¥ +. + Bes 1)| 


CA Po—P, 


Substituting in this the value of Y and using the BET equation (63) to change 
the independent variable from N to P, we get 


pets 


Dp Sw= Nx. (73) 





Ja tr CP? (48 — 2)CP 
CA Po— Py, \4 (1 —P)\1—P+CP)P? (1—P)\(i—P+ CP) 


— (28 — 1) log | B4/ CP? 4: (48 — 2)CP 
(1 -P)1-—P+CP)P (1—P)(l1—P+CP) 





P 





\}P=Po 
ato Se lead a. + B(2e6 — »|} : (74) 


P)(1 — P+ CP) Jr-P, 
or reverting to our former notation where 
Waco (ome, by. (48 — 2)CP 
ti ~ PL ~F 4 OPyp tl ~ Ph CP) 


we have 


l |, a 


\ 


— BCP» 

; —_—______—_——___—. + B(26 - 1 
Vita ae ee Ce )} 
Vy_t+ 


+ B(26 — 1) 





: BCP, 
(1— P,)(1 — Pz+CP,) 


| Re 


/ 


The value of D, / B’kT given by Equation (75) has been plotted in Fig. 5 as 


, ; ps aa to s : 
a function of the mean relative humidity te, In this calculation P; has 


been taken as zero. The values of B and C have again been taken from the 
data on spruce wood (1) (C = 11.80 and B = 6.17). 


There are not at present sufficient measurements to test this equation 
quantitatively but the measurements that are available do indicate that the 





BABBITT: THE DIFFERENTIAL EQUATIONS OF DIFFUSION 471 


experimental relation between D, and P is represented qualitatively by Equa- 
tion (75). All measured results for wood seem to indicate a diffusion coefficient 
very small at low humidities and increasing very rapidly at humidities ap- 
proaching saturation (2, 8, 11). A curve of similar shape has been obtained by 
King (7) from his measurements on the diffusion of water vapor through horn 
keratin. In fact, all available measurements indicate that this shape of curve 
is general for diffusion whenever the adsorption can be represented by the 
sigmoid isotherm peculiar to BET adsorption. 


———L_ 4 —t — 
O2 o3 O4 Qs 
MEAN RELATIVE HUMIDITY 


Fic. 5. Variation of the coefficient of the integrated equation when the adsorption obeys the 
BET equation. 


The part of the theoretical curves of Figs. 4 and 5 that has not been observed 


experimentally is the reversed curvature and sharp decrease to zero at low 
relative humidities. In the experiments, however, few measurements have 
been made at the lowest humidities and it might well be that this reversed 
curvature has been overlooked. More measurements at low humidities will be 
necessary to delineate the experimental curve in this region. The theoretical 
curves do show, however, the large increase of diffusion at high humidities 
that is so striking a characteristic of diffusion through-adsorbing solids. 


It should be pointed out that these diffusion equations have been derived 
on the assumption that the adsorption can be represented by the BET equation. 
It has been shown by numerous workers that the BET equation gives satis- 
factory values for the amount of gas adsorbed only at relative humidities 
below 50% and therefore we should not expect results derived from this theory 
to hold outside of that range. Although it does not give quantitative agreement 
at high humidities, the BET equation does indicate, however, the large increase 
of adsorption that is so characteristic of the adsorption of water vapor as 
saturation is approached and, in consequence, the equations can be considered 
as giving a qualitative indication of the behavior at high humidities. 


Flow Through a Medium that is Both Porous and Adsorbing 


There remains for consideration the problem of a medium that is porous in 
the sense that it is permeable to a nonadsorbed gas and that is at the same 
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time an adsorbent for the gas whose flow is being measured. It is obvious that 
with such a medium there will be a movement of the gas both through the 
pores and through the material itself. It is also obvious that through the pores 
the quantity of flow will depend on whether we have hydrodynamical flow 
with an absolute pressure difference or diffusive flow with a difference only of 
partial pressure, whereas through the medium proper, as shown above, the 
flow will depend only on the gradient of ¢ and will be disturbed by the presence 
of the second gas only indirectly through the altered pressure gradient in the 


pores. 


Although at first sight one might suppose that such a flow could be repre- 
sented by the sum of the mass transfer through the pores and through the 
medium, this is not so since in general there will be throughout the medium a 
free exchange of gas from the adsorbed phase within the medium to the vapor 
phase within the pores. In a small element of internal surface there will be a 
balance between the mass flowing into the element, the mass flowing out, the 
mass that evaporates into the pores and the mass that condenses on the surface 
from the vapor in the pores. Similarly in a small element of volume in the pores, 
the adsorption and desorption of the gas by the medium must be taken into 
consideration in addition to the flow. 


Mathematically, this amounts to the fact that it is not valid to apply the 
equation of continuity separately to the diffusion through the medium and 
to the flow in the pores, but instead the equation of continuity must only be 
applied to the total mass transferred, that is, to the sum of the mass trans- 
ferred by diffusion through the medium and by flow through the pores. Thus 
although it is simple to set up the fundamental dynamical equations for the 
flow it is not possible, in general, to obtain from them the differential equations 
since the equation of continuity could only be written if we knew precisely the 
size and distribution of pores in the medium and could calculate at each point 
the distribution of gas between gaseous and adsorbed phase. We shall not at 
present attempt to do this but shall rest content with having pointed out the 
difficulties. 

Conclusion 


We have in this paper derived from first principles the differential equations 
governing the diffusion of gases and vapors and have endeavored to point out 
the modifications that must be made to these equations when the diffusion is 
through a solid medium. In particular we have shown that, where there is an 
interaction between the diffusing gas and the solid, as in adsorption, the use 
of Fick's law as the fundamental equation of diffusion is not valid. An attempt 
has also been made to distinguish between hydrodynamic flow and diffusion 
and to point out the assumptions and features peculiar to diffusion. 


In conclusion, a general comparison of the derivation of the differential 
equations of four equalization processes is given in Table I. In this table the 
equivalent fundamental quantities are set down and from these quantities 
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step by step the analogous equations are derived. Such a comparison is useful 
to show where the processes differ and whence the variations in the form of the 
final equation. 


There are one or two points that are brought out by this table that are 
worth emphasizing. In the first place it is essential to note that the pressure 
(and its analogous quantity @) is the potential function whose space derivative 
gives the force that is balanced by the resistive forces. This fundamental fact 
0c 
Ox 
as the fundamental equation of motion. As brought out in this paper, only 
when ¢ is directly proportional to the pressure is such an approach correct. 


has been obscured by the tendency to use Fick’s law in the form P= — D, 


TABLE I 


Conduction of | Flow according to } Diffusion of | Diffusion through 
| heat Darcy's equation | gases | absorbing media 
Fundamental |Energy Mass Mass | Mass 
quantities |Quantityof heat =|Density = mass perce. |Density = mass}|Number of adsorbed 
| energy perce. =q| =p per cc. = p | molecules per cc. = N 
Velocity of flow of Velocity of flow = 7 |Velocity of diffu-| Velocity of diffusion = u 
energy =v |} sion = 4 | 
Flux of heat = f =|Flux = pv Flux = rate of [Flux = rate of flow of 
| qv | flow of gas = pu | adsorbed gas = Nu 
Temperature = 6 |Pressure = p Partial pressure {Spreading pressure = ¢ 
| =p | 


| 
| 
| 


a i eee 5 | lo “ ae 
? “ nts 5 servati | : : ¢ | »ynservation o . : 
Fundamental | Conservation of | Conservation of mass | COmServation of | Conservation of mass 
law | energy mass 


Statement of | Paani 00 | gs pi(u;—u,) = i dg, 
fundamental i ax | : u ox 1 apr * Ox 


a iwhere k = permeability Cp2 ax where C = coefficient of 
- | age eel 90 ww = viscosity ; .. | resistance 
> ae’ | ‘ where C = coeffi- 
I wal il 3 att cient of resistance 
where k = thermal) 
conductivity 


Equation of | a(f) 2 Oe O(pu) La OP, | 9(p1%1) ke a(Nu) 
continuity | “gx I ax at Serr j Ox 
jor |where f = porosity 
O(qv) _ 4q | 


ax at 





Equation of | qQ = pcd, p = pop™ pi = Rif ¢ _ F(N) 

state where p = density|[sothermal flow: m (flow assumed The relation between ¢ 
C =. OG eg isothermal) and N is dependent on 

fic |Adiabatic flow: m = v the nature of the adsorp- 

heat | tion and is related to the 

adsorption isotherm 


j 
| 
j 
| 
| 
| 


Differential oq _ ,, 90 kpo apm+i ees en... 3-8 (w a) 
cquaseu? at ax? = u(m+1) dx? ; at , C2 Be 
00 =k «0 
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_ 86 
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The second point that should be noted is that the dynamical equaticn has 
two essentially different forms; in the conduction of heat and in the inter- 
diffusion of two gases, the flux is proportional to the potential gradient, whereas 
in flow according to Darcy’s equation and in the diffusion in an adsorbing 
medium it is the velocity of flow. Since in all four processes the equation of 
continuity pertains to the flux, this difference in the dynamical equation leads 
to a difference in the final equations. 


It would scarcely be necessary to mention the effect of the equations of state 
were it not for the confusion that has arisen from the indiscriminate use of 
Fick’s law. The classical example of the influence of the equation of state, 
although never expressed in that way, is the difference between the thermal 
conductivity and the thermal diffusivity as shown in the different forms of the 
final differential equation for the flow of heat. If we look at the analogous 
equations in the interdiffusion of gases, we find that the diffusion coefficient 
D is the counterpart of the thermal diffusivity x whereas Dk, is analogous to 
the conductivity. It is necessary, therefore, to keep these relations constantly 
in mind so that the pressure and the density are not indiscriminately used in 
the equations of diffusion. 


Finally it is essential to point out that this approach to the diffusion through 
adsorbing media is tentative, pending greater confirmation than it has here 
been possible to obtain; the assumption that C is independent of the number 
of molecules adsorbed may not be valid and it may be necessary to take C as 
some function of NV. There is no doubt, however, that the true differential 
equation for the diffusion of a gas in such a medium can result only from a 
combination of the equations of continuity and of state with an expression 
defining the fundamental dynamical law of the system. 
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